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A review of the present state of investigations of the pseudospin-electron 
model (PEM), which is used in the theory of strongly correlated elec- 
tron systems, is given. The model is used to describe the systems with 
the locally anharmonic elements of structure represented in the model by 
pseudospins. The consideration is based on the dynamical mean field the- 
ory approach and the generalized random phase approximation. Electron 
spectrum and thermodynamics of the model are investigated; the cases 
of the simplified model, the model with strong interaction and the two- 
sublattice model are studied more in detail. The phase transitions into 
other uniform or modulated states as well as superconducting phases are 
described; the criteria of their realization are established. Based on this, 
the description of structural and dielectric (ferroelectric type) instabili- 
ties, phase separation and bistability phenomena is given. A comparison 
is made with the thermodynamics of the Falicov-Kimball model (which 
can be considered as a particular ceise of PEM). The possibility of ap- 
plying the PEM to the analysis of thermodynamics of the real HTSC 
systems is discussed. Attention is paid to the unsolved problems in the 
study of PEM. 
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1. Introduction 

Much attention is paid in recent years to the investigation of systems with 
strong electron correlations, such as crystalline compounds with the transi- 
tion and rare-earth ions (transition metals, transition metal oxides, mixed 
valent compounds, heavy fermion systems, high temperature superconduc- 
tors, etc.). Their specific electronic, magnetic and conducting properties as 
well as the presence of a variety of phase transitions and the phenomena 
that are connected with this property are caused to a great extent by the 
splitting and reconstruction of energy spectrum due to correlation effects. 
The theory of such systems is based on the Hubbard model and on its gen- 
eralizations, where the crucial idea is the one concerning the decisive role 
of the strong short-range interaction of particles (electrons). At the pres- 
ence of additional (i.e., vibrational) degrees of freedom, one can mention, 
among others, the pseudospin-electron model (PEM). The model appeared 
recently in connection with the investigation of the high- Tc superconduc- 
tors. It was introduced to describe the contribution of locally-anharmonic 
elements of the crystal structure to their electronic properties. 

Electron system in PEM is described by the Hubbard Hamiltonian while 
the anharmonic vibrational modes are treated using the pseudospin formal- 
ism. The model Hamiltonian is as follows 

H = J2[Un^,^n,,_l + {gS^-^l)in,,^+n,^l)-hS:-nS^] + Y, Uja+a^^j.,,. (1) 

Here at^^, o.f^ are electron annihilation and creation operators, rii^a is an 
electron occupation number; besides the electron correlation (J7-term), the 
single-site part includes the interaction with pseudospin (g-term) and the 
energy of the tunnelling- like splitting of vibrational levels (fi-term); the 
field h describes the asymmetry of local potential. The electron transfer 
(t-term) is included as well. 
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The pseudospin-electron Hamiltonian (in form ([T]) with inclusion only 
of the ^-interaction) was used by Miiller with the aim of describing the 
anharmonic vibrations in the oxygen subsystem of the high-Tc supercon- 
ducting crystals of the YBaCuO typeffl The YBa2Cu307_5 crystal is a 
typical and most studied example of such objects. The unit cell contains, 
besides two superconducting planes, the chain (at the (5 <C 1 composition) 
elements Cui-Oi, connected by CU1-O4-CU2 bridges with CU2-O2 plains 
through the apical oxygen ions O4. The vibrations of these ions along the 
c-axis (perpendicularly to the plains) exhibit a strong anharmonicity. Much 
evid ence exists in support of this concept. One can mention the EXAFS 
data,'212l Raman s catteri ng and dielectric measurements, SEI^EIHini local po- 
laron phenomena! ^ ^ 1 ^ -^ J bistabilities in the normal phase regioJI^ as well 
as neutron scatteri ng i nvestigations or the results of the first principle 
LAPW calculations .IJ^ Despite a certain ambiguity in the data, the conclu- 
sions were made about the existence of two different equilibrium positions 
of the O4 ion. The local double-well picture as an approximate simplified 
model was supported by the oxygen O4 vacancy effect on the positions of 
apical ions observed in Ref. T5]. 

Moreover, a connection between positions of O4 ions and electron states 
in CU2-O2 plains plays an important role in YBa2Cu307_5 crystals. The 
data given in Ref. 16J point to the existence of a significant correlation 
between the occupancy of electron states of the Cu2 ion and the R04-CU2 
distance as well as to the decrease of this distance at the transition from 
the metallic orthorombic phase to the semiconducting one (that takes place 
a.t S > S* = 0.55). These and other similar facts suggest the presence of a 
large electron-vibrational coupling. In the pseudospin representation, when 
the pseudospin variable S^ — ±1/2 defines the positions of O4 ions, it is 
described by the gS^rii^ term. 

Consideration based on the PEM was applied to the HTSC systems 
starting from Refs. [I|17l - f20] . Hamiltonian similar to ([T]) was used by Hirsch 
and Tang in the study of electron states in the framework of cluster cal- 
culations. In the context of the idea concerning t he effect of a nharmonic- 
ity on the superconducting transition temperature,^l^^I^^I^^2I] the possible 
connection between supercon duc ting pairing and the lattice anharmonicity 
was considered by Frick et al.^^ (the quantum Monte-Carlo calculations). 
In what follows, the investigations of the PEM were devoted to the a nalysis 
of the electron spectrum,!^ the pseudospin and collective dynamics ] I I 
the charge and pseudospin pair correlations and the behaviour of dielectric 
susceptibility.'^ 
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In subsequent investigations, the main attention was paid to the ther- 
modynamics of the model in special cases and simplifications: (i) a model 
with the infinitely large correlation (U — > oo) when the double occupa- 
tion of electron states on the site is excluded; (ii) simplified PEM with 
U = and = 0; (iii) simplified model {U = 0) with the tunnelling- 
like dynamics [il ^ 0); (iv) two-sublattice PEM for the layered structures 
of the YBaCuO- type; (v) the cluster PEM of ferroelectric-superconductor 
heterostructures! 

There was performed a study of phase transitions between the states 
with different electron concentrations and with different orientations of 
pseudospins (in the regime of the fixed chemical potential, fi = const), 
and the p hase se paration effects at a given concentration of electrons 
(n = const) !2213Ql The possibility of the appearance of a doubly modulated 
(so-called chessboard) phase o r an inco mmensurate phase (in the case of 
weak coupling) was 

establishedpmiSS] 

the superconducting instability in 

OA 

PEM was analysed;^ In the case of two-sublattice PEM the structural in- 
stabilities of the ferroe lectric type as well as the bistability phenomena were 
revealed and analysed P^^^^^^ZEH xhe study of the PEM thermodynamics 
was performed mainly within the generalized random phase approximation 
(GRPA)!321 A method of dynamical mean field theory (DMFT) was used in 
the case of simplified PEM,^ when the analytic formulation of the theory 
is possible. 

As was shown, the PEM also possesses an interesting collective dynam- 
ics of pseudospins. The corresponding spectrum change s its form depending 
on the electron concentration and temperatureP^ES ^^le spectrum is dif- 
ferent at high or small values of g and its shape also depends on the h and 
parameters. In Refs. ^40ri42) . the contributions into Raman scattering 
intensity, connected with the mentioned collective pseudospin excitations 
(that correspond to the phonon-like vibrations of anharmonic subsystems in 
the YBaCuO structures) and electron intraband and interband transitions, 
were considered. 

The PEM is closely related to the Falicov-Kimball (FK) model inten- 
sively studied in recent years (see, for example, Ref. 03]), in which the 
interaction between the localized and itinerant particles (electrons) is re- 
sponsible for the similar phase transitions (between states with different 
concentrations of particles and with/ without spatial modulation). The sim- 
plified version of PEM corresponds to the FK model in case there is no 
tunnelling-like splitting in the PEM and when the localized and the mov- 
ing particles in FK model have different chemical potentials (one can pass 
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to FK model putting — Wi — 1/2 and h = jl, where Wi and fi are occupa- 
tion number and chemical potential of localized particles, respectively). The 
regimes of thermodynamic averaging are usually different for both models 
(a fixed concentration of the localized particles for FK model and a given 
value of the field h for the PEM). 

It should be mentioned that for the recent few years the PEM has found 
application in describing the cha rge transfer in molecular and crystalline 
systems with hydrogen bondsP^The model is also very promising in inves- 
tigating the thermodynamics of processes connected with ionic intercalation 
in the layered structures (see Ref. [IS]), where different separated positions 
exist in the unit cell for the intercalate ion, and the hopping between them 
is possible (intercalation of Li"*" ions in the Ti02 matrix provides an exam- 
ple of such a situation^Sl). A specific version of PEM was recently usecP^ in 
modelling the electronic properties and the field effect in the Cu02/SrTi03 
interfaces in HTSC/STO heterostructures. 

This paper presents a review of the main results concerning the thermo- 
dynamics and energy spectrum of the PEM in the above mentioned cases 
and approximations. The case of the simplified PEM {U — 0; Q, — Q or 
0) is considered more in detail. Dynamic properties of PEM are not 
considered here. An attention is paid to the possible application of the 
PEM to the description of inhomogencous states, structural instabilities 
and bistability phenomena as well as transitions into the phases with the 
charge modulation in the high T,. superconductors and other systems to 
which the model can be applied. 

2. Thermodynamics of Simplified PEM in Dynamical Mean 
Field Theory 

The dyn am ical mean field theory approach proposed by Metzner and 
VoUhardiPSl for the Hubbard model (see also Ref. [JS] and references 
therein) is a nonperturbative scheme which is exact in the limit of the 
infinite space dimension [d — > oo). The method is very successful in con- 
sidering the systems with strong electron correlations and is used with ad- 
vantage in solving a variety of problems and models. Within the framework 
of DMFT, investigations of the single-particle spectrum and the thermo- 
dynamics of the simplified PEM were performed for a strong coupl ing case 
{g ^ W, where W is the half- width of the initial electron band).'22l The 
Fourier-transform Ga {ujn , k) of the electron Green's function 
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(where T denotes the r- ordering procedure) with the scattering matrix 

a{(3) =Texpi^- J^tJ^ U, a+ {T)a,, (r) | (3) 



and the averaging with the single-site part iJp ~^Hi oi the Hamiltonian 

i 

^ , is expressed as a series in terms of the electron hopping parameter Uj . 
The Larkin's equation 

Gr,(r - r') = Sr,.(r - r') + E^r - t")U^GI,^{t" ~ r') (4) 

separates the total irreducible (with respect to tij) part S"^; formally 

G.(^n,k)= „_ \ . (5) 

In the case of high dimensions (d —^ oo), when the hopping integral is scaled 

{tij tij/Vd), only single-site contributions survive in the expression for 
l50l 

E'ljir - t') ^ SijE^{T ~ t'); E^{ujn,'k) ^ E^{ujn). (6) 

Such a site-diagonal function, as it was shown by Brandt and MielschjSH 
can be calculated by mapping the infinite-dimensional lattice problem on 
the atomic model 

e-PH ^ ^ ^-mo (7) 

xTexp|-y dry dr' Mr - T')at {T)a„iT') 

with auxiliary Kadanoff-Baym field Jo-(t — t'|^ which should be selfcon- 
sistently determined from the condition that the same function S^- defines 
the Green's functions for lattice ^ and atomic limit 

G(^\u„)^ _ y (8) 

"Dynamical" mean field Ja{T — r') (so-called coherent potential) de- 
scribes the hopping (transfer) of electron from atom into environment at 
the moment t, and propagation in environment without stray into atom 
until moment r'. The connection between this "dynamical" mean field of 
atomic problem and Green's function of the lattic e ca n be obtained using 
standard coherent potential approximation (CPA)E21 

J,(^„) =S;i(c^„)-G;i(a;„), (9) 
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where 

gWk) = g.K)= r°^dt _^f\ (10) 

is a single-site Green's function both for atomic limit and lattice. Here 
summation over wave vector was changed by the integration with the 
density of states (DOS) p(t) (the Gaussian one for the hypercubic lat- 
tice p{e) — -j^T^e"^ and semi-elliptic DOS for the Bethe lattice 

p{e) = ^^W^-e^ see Ref. [49]). 

In order to find expression for Green's function in the atomic limit, one 
can use the fact that the statistical operator of the single-site problem ^ 
can be expressed in the forrrPSI 

e-,9//eff = p+^-0H+ + p-e-0H- (11) 

because the atomic space of states splits into two independent subspaces. 

As a result, the single-electron Green's function is a sum of Green's 
functions in subspaces and is equal to 

Gi"^ K) = ■ ^ ^^^l , g + . ^ ^^'l ■ (12) 

lUJn + P - Ja[i^n) - i ZCJ„ + /i - J^r (tJ„ ) + | 

Here (. . .) is the statistical averaging with the effective Hamiltonian ([7|). 
Partition functions in subspaces are 

Z± =Spe-'3^± ^ei'^-Qi (13) 



e±^n 1+^^^^"^'^ n 1- 



Pseudospin mean value is determined by the equation 

^^'^ ^ I zl + f = ^tanhi (/?/.- (g+[(5^)] - Q_[(5^)])) . (14) 
Electron concentration mean value is determined by 

and the functional of the grand canonical potential can be derived in the 
standard way for DMFT 



N 



'^l-) - ^ E 1 1^ G-"^ ('^^ ) - 5] In (c.„ , k) I , 



(16) 
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where 



(a) 



Z-) 



(17) 



is a thermodynamic potential for atomic problem. 

The solution of the above given set of equations and the calculation of 
thermodynamic potential were performed for the case of semi-elliptic DOS. 
The field Ja{LOn) is determined by the simple cubic equation 

(P+) , (p-> 



M/2 



lUJn 



JaM + f 



(18) 



/i - Ja(t^,0 - f ILUn+l-l- 

The solutions with 3mJo-(cL;) > are considered; the condition 3mJcr(cj) — > 
determines the band boundaries. Their dependence on coupling constant 
at the fixed value of (S^) is shown in Fig.[TJ It can be seen that there exists 



Im J(u)=0 




Im J(a)>0 






Im J(u)=0 


Im J(a)^0 





1.5 g 2.0 



Fig. 1. Electron bands boundaries (semi-elliptie DOS, W = 0.4, (S^) = 0.2). 



a critical value oi g ^ Q.bW when a gap in spectrum appears and we have 
an analogue of the Mott transition. In the case when the single-electron 
Green's function is calculated in Hubbard-I approximation (the scattering 
processes via coherent potential are not taken into account, Ja{u)) = and 
S^(w„) = (P+)/(«a;„ + ^l- g/2) + (P-)/(iu;„ + + g/2)), the electron 
subbands are always split and the gap in spectrum exists at any values 
of g (see below). From this point of view, the Hubbard-I approximation is 
insufficient; even in the case of strong coupling {g 3> W) it only qu alitatively 
describes the dependence of the subband half- widths on {S^)\^ 
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The expressions presented above allow us to investigate in the DMFT 
approach the thermodynamics of the simplified PEM. It was done in 
Ref. 26J in the /i = const and n = const regimes. 

In the first case, the thermodynamically stable states are determined 
from the minimum of the thermodynamic potential (fT6|l . Analysis of so- 
lutions of CPA equations for Ja{uin) together with the equation for 
{S^) shows that in this regime the first order phase transitions with the 
jumps of the pseudospin mean value and electron concentration can take 
place. Such transitions are realized when the ^ and h values correspond to 
the split subbands in an electron spectrum (see the phase diagram (/i — h) 
at T = in Fig. [2). The field dependencies of (5^) and grand canonical 











M : 


S'=- 


1/2 


n,=2 




n= 


z 




0.5 - 




0<7i<2 










-0.5 - 


0<7i<2 










T1 = 






71 = 




S'= + 1/2 



(a) 




Fig. 2. Phase diagram (^t — h). Dashed and thin soUd hnes surround regions with 
5^ = it^, respectively. The lines of the first order phase transition are shown in bold, 
a) g = 1, = 0.2; b) g = 1, = 0.7. 



potential in the region of the phase transition point are shown in Fig. [3l 
Since the band structure is determined by the pseudospin mean value, the 
change of the latter is accompanied by the corresponding reconstruction 
of the electron spectrum. With the temperature increase the region of the 
phase coexistence narrows. The corresponding phase diagram (Tc — K) is 
shown in Fig. [H One can see that with respect to the Ising model the phase 
coexistence curve is shifted in the field and deviates from the vertical line. 
Hence, the possibility of the first order phase transition with the tempera- 
ture change exists in the pseudospin-electron model for the narrow range 
of h values. 
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(a) 



(b) 



Fig. 3. Field dependencies of {Sz) (a) and grand canonical potential (b) for ij, = const 
regime when chemical potential is placed in the lower subband /i = —0.37 {W = 0.2, 

9 = 1, r = o). 
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Fig. 4. Phase diagram (Tc — h): solid and dashed lines indicate the first order phase 
transition line and boundaries of the phase stability region, respectively (g = 1, W = 0.2, 
At = -0.5) 



When the electron concentration is fixed (regime n — const), the first 
order phase transition transforms into the phase separation. The regions 
appear where the derivative d^ijdn is negative. Figs. [5] and [6l The corre- 
sponding phase diagram (T — n) is built (Fig. [7l see also Ref. [26] ) with the 
use of the "Maxwell rule" which follows in this case from the replacement of 
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the original free energy in its concavity region by the tangent hne. The dia- 
gram describes the separation to the states with the large and small electron 
concentrations (and with the {S^} w —1/2 and (S^) w +1/2 pseudospin 
averages at low temperatures, respectively). 

1 ,0 n 



0.0 - 



-0.5 




-1.0 I 

0.0 0.5 1.0 1.5 2.0 

n 

Fig. 5. Dependence of the chemical potential fi and electron bands boundaries (dashed 
lines) on the electron concentration n (T = 0.001, g = 1, W = 0.2, h = 0.1). 

It should be noted that the problem of phase separation in strongly 
correlated systems is not new (see R ef. 1 53j and references therein). It was 
shown for Hubbard and t— J modelJ^SI tJiat for some values of the param- 
eters, the system separates into hole-rich and hole-poor regions with para- 
magnetic and antiferromagnetic orders, respectively. In our case of PEM 
without electron correlations, the system separates into regions with elec- 
tron spectrum that contains both wide empty electron band and occupied 
localized states (at n ~ 0) and partially filled wide electron band and 
empty localized states (at n ~ 1), see Fig. [S] their weights are determined 
by the electron concentration. Localized states of such a type (polarons) 
result from the strong electron-pseudospin coupling (strong interaction of 
electrons with the out of plane apical oxygen vibrations) in the case of 
YBaCuO-type structures, and it can be supposed that the hopping be- 
tween such polarons manifests itself in the carrier relaxation.!^ 

For the first time the possibility of phase separation in PEM was men- 
tioned in Ref. [56] where it was considered within GRPA in the limit of 
strong correlation U ^ oo (see below. Sec. |4]). Here it is obtained for the 
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Fig. 6. Dependence of the chemical potential ^ (a) and the deviation of free energy from 
linear dependence AF = F{n) — ^F{2) — (l — f^) ^(0) (b) on the electron concentration 
n for different temperatures T {g = 1, W = 0.2, h = 0.1). 




Fig. 7. Phase diagram (T — n) for phase separated state. Solid line: binodal, dashed 
line: spinodal {g = l,W = 0.2, h = 0.1). 



opposite case of [/ = 0. As a whole, such a picture of phase transitions into 
a new uniform phase is in agreement with the known resuhs for the FK 
model in the case of strong coupling (see Ref. [43]). But, as is evident from 
the phase diagram obtained in Ref. [S^ for this model, in the region of 
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large but finite values of g the phase with double modulation can appear. 
In order to detect instabilities associated with the wave vector k 7^ one 
should calculate the susceptibility functions and analyse their temperature 
and k- dependencies. Such an investigation was performed for PEM in the 
framework of GRPA (see below, Sees. [3X2l and [3X2]) . 

3. Simplified PEM in Generalized Random Phase 
Approximation 

A more complete investigation of the PEM was performed in the framework 
of the generalized random phase approximation (GRPA). Such an approach 
was formulated by Izyumov and Letfulov for the calculation of the pair 
correlation functions and magnetic susceptibility of the Hubbard and t — J 
models. It is based on the expansions in terms of electron transfer and 
consists in the summation of the diagrams having a structure of sequences 
of the electron loops (created by the electron Green's functions) j oined b y 
vertices of various types appearing due to short-range interactions.^^^^SSl 

Thermodynamics of PEM was studied in the GRPA ior U = 0, fl = 
(simphfied model) in the cases of strong {g W) and weak {g < W) cou- 
pling as well as in the limit of the infinitely large on-site electron repulsion 
(U — > 00). Let us first consider the results obtained for simplified model. 

3.1. Strong Coupling Case; U — 0, D, — 
3.1.1. Thermodynamics of the Uniform State 

First we consider the case of strong coupling. Here the single-site states 
can be used as the basic ones and the formalism of electron annihilation 
(creation) operators Oia- — hic^P^ , — birrPf {Pf" — 1/2 ± S^) acting 
at a site with certain pseudospin orientation was introduced.^^Hl Using this 
representation we can write the model Hamiltonian in the form 



Here, e = — /i + g/2 and e = — /i — .g/2 are the energies of single-site states. 

Expansion of the calculated quantities in terms of electron transfer leads 
to the infinite series of terms containing the averages of the T-products of 



H ^ Hq + Hint , 





(19) 
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the Oicr, o,i„ operators. The evaluation of such averages is made using the 
corresponding Wick's theorem. The results are expressed in terms of 
the products of nonperturbed Green's functions and averages of a certain 
number of the projection op erat ors Pf^ which are calculated by means of 
the semi-invariant expansion.!^ 

Nonperturbed electron Green's function is equal to 



P7 



ILO,-, — £ ILU„ — £ 



(20) 



In the diagrammatic representation it has the meaning of the simplest ir- 
reducible Larkin part in the series for the free single-electron Green's func- 
tion Gk(wn). In the Hubbard-I type approximation (see previous section) 
Gk(wTi) can be written as a sum of the following chain diagrams 




or 



Gk(c^«) = [g \w,i) - tk] 

and its poles (after analytic continuation zcij„ 
electron spectrum 



. , (21) 
(22) 

iS) determine the 



£iji{t\i 



2 



1 



^i±-Jg-^ +At^{S^)g + i 



k • 



(23) 



The electron subbands are divided by a gap which tends to zero only at 

There are used here (and below) the following diagrammatic notations: 
Q—S^, — Cl^^^^O — ^9i{^n), wavy line is the electron intersite hopping Uj. 
Semi-invariants are represented by ovals and contain the (5-symbols on the 
site indices. 

In the adopted approximation the diagrammatic series for the pseu- 
dospin mean value can be presented in the form 




(24) 



In the spirit of the traditional mean field approactP^ the renormalization 
of the basic semi-invariant by the insertion of independent loop fragments 
is taken into account in (l24ll. 
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The analytical expression for the loop is as follows: 




-y— 



n,k 



P7 



Pr 



l3{aiP+ +a2P-). 



(25) 



Similarly, the diagrammatic series for the electron concentration mean value 
is as follows: 



=(Q) 




where Q —Ui 



(26) 



£" = (£,£), 



ptHP^.pn- 

The grand canonical potential $ and pair correlation functions {{S^Sj), 
{S^rij), {niUj)) are calculated according to self-consistent scheme of the 
GRPA: in sequences of loop diagrams in the expressions for $ and correla- 
tors the connections between any two loops by more than one semi-invariant 
are omitted. We have, respectively 




(27) 



(28) 
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(29) 



The first term in equation ([28|) takes into account a direct action of the 
internal effective self-consistent field on pseudospins: 



(30) 



leading to the renormalization of the second-order semi-invariant due to 
the inclusion of "single-tail" loop- like parts. Second term in equation ((28l) 
describes an interaction between pseudospins which is mediated by electron 
hopping. From p4|) and (|25p . the equation for the pseudospin mean value 
follows: 




(5^) = ltanh|^(/. 



a2- ai) 



In 



1 



-/3£ 



1 + c-/3g 



(31) 



where 



a2 — ai 



L X ^ i. 



Ke//(ik)) -n(e,(tk))] . (32) 



The grand canonical potential in the considered approximation has the 
form: 

(cosh |£jfa))(cosh f £//(tk)) 



A$ = $ - $ 



t=0 



Ni3 



k 



+ (5^)(«2 



1 

7' 



■ In cosh 



/3 



1 



In- 



r/3 

^ In cosh { —ih 
13 \2^ 



(cosh |e)(cosh fe) 
ai) 4- In 



2 o y i^v^vjoii 2 

1 + e~P^ 



Ol2 



1 e-^^^ 



/i + ln- 



1 



-/3e 



1 + , 

With respect to the initial GRPA schemepSEH action of the internal 
effective self-consistent field on pseudospins is taken into account by includ- 
ing the mean field type contributions into the expressions for all thermody- 
namic quantities. The electron concentration and pseudospin mean values 
as well as correlation functions are calculated co nsisten tly with the thermo- 
dynamics functions. It can be checked explicitlj^^nEI] using the relations 
d$ d$ dS*^ 



d(-M) 



d{-h) 



di-ph) 



(33) 
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At high temperatures, equation ([?T|) possesses only a uniform solution 
(Sf) = (S^). However, there exists a possibility of phase transition be- 
tween different uniform phases with the different pseudospin mean values. 
For the first time the possibility of such a transition, which leads to the 
structural (dielectric) instability, was considered for the PEM in the limit 
of the strong electron correlation {U oo) in Refs. |25|26j (this issue is 
analysed below, see Sec. S]). A relatively complete description of this tran- 
sition was also given in Ref. [SD] for the PEM with a direct interaction 
between pseudospins (in the tij = limit). 

Within the GRPA scheme presented here, the uniform-uniform phase 
transition in the simplified PEM was analysed in Ref. [3D] . The solutions of 
the set of equations (j33p which correspond to the absolute minimum values 
of A$ were determined and analysed. The calculations were numerically 
performed for the square lattice with the nearest-neighbour hopping (d — 2 
DOS with the bandwidth 2W). The obtained picture of the phase transition 
in the /i = const and n = const regimes is very similar to the case of DMFT 
approach. When the chemical potential is fixed, there exist jumps of the 
pseudospin mean value and electron concentration on the phase transition 
line (Fig. [8]). The {h — fjt) and (Tc — h) phase diagrams are nearly the 
same by their shape as in the DMFT case despite the fact that electron 
energy spectrum in GRPA is split at any values of g, see Fig. [9l [TOl The 
same conclusion can be made when we compare the (T — n) phase diagrams 
obtained in the n = const regime; the corresponding phase separation areas 
are shown in Fig. [TlJ The difference in positions of critical points in the (/i, 
T) plane calculated in the DMFT and GRPA approaches is also small (see 
Figs. H and [lOb). 

In figures discussed above, the case is presented when the chemical po- 
tential is placed in the lower energy subband. If ^ is placed in the upper 
subband, the results are transformed according to the internal symmetry 
of the Hamiltonian and the following replacements should be made 

H^-H, h^2g-h, n^2-n, 5^ ^ -5^ (34) 

Sur nming up, we can conclude that the generalization of GRPA 
scheme,!^^^^ which takes into account the mean-field loop-like contribu- 
tions to the semi-invariant averages, makes it possible to calculate the ther- 
modynamic functions and to investigate the first order phase transitions 
between different uniform phases. At the transitions, there always remains 
a gap (in the strong coupling case, g W) in the electron spectrum. With 
the change of the mean value of the pseudospin a reconstruction of the elec- 
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Fig. 9. Electron bands boundaries (to = 0.4, (5^) = 0.2). 



tron spectrum takes place, at which the widths of the electron subbands 
change which results in the jump-like change of electron concentration (at 
the given chemical potential fi, it corresponds to a charge transfer from /to 
the electron reservoir). 

The phase coexistence curve in the (h, T) plane is tilted from the vertical 
line; therefore, there is a possibility of the first order phase transition with 
the temperature change (in the narrow interval of the field h values). It 



2, 2008 3:3 



WSPC/Trim Size: 9in x Gin for Review Volume 



stasyuk'new 



Phase Transitions in the Pseudospin-Electron Model 19 




(a) (b) 

Fig. 10. (a): {h — fi) phase diagram (T = 0, to = 0.2, g = 1). (b): (Tc — h) phase diagram 
(g = 1, to = 0.2, fi = -0.5). 




Fig. 11. (T - n) phase diagram {h = 0.1, tfe^o = 0.2, g = 1). 

should be noted that the existence of the shifted and tihed coexistence curve 
(as the result of the local pseudospin-electron interaction) was obtained for 
the first time in Ref. [60] for a PEM with the direct interaction between 
pseudospins. In the case considered here, such a direct interaction is not 
included, but due to electron transfer there appears an indirect one (the 



February 2, 2008 3:3 



WSPC/Trim Size: 9in x Gin for Review Volume 



stasyuk'new 



20 /. Stasyuk 

latter is formed by the loop- like contributions 1 11 | q). The relative role of 
the direct and indirect interactions between pseudospins is analysed more 
in detail in Sec. [SJ on the example of the two-sublattice PEM. 

Herein below, considering the (S^S^) pair correlation function, we shall 
look at the possibility of the spatially modulated charge (charge density 
wave, CDW) and pseudospin orderings in the PEM with U — 0, — in 
the strong coupling case. 

3.1.2. The Chess-Board Phase 

The analysis of the {S^S'^)q correlator temperature behaviour shows that 
for certain values of model parameters, the high tem perature phase may be 
unstable with respect to fluctuations with q ^ 0.^^ Solution of Eq. ((28|l for 
pseudospin correlator has the form 

1/4- (5^)2 



1 



n 



where | n | q characterizes an interaction between pseudospins via electron 
subsystem: 



n 



1, , (36) 

a,/3 



and its singularities 

'n 



,(7-(^^>') = 



give the instability points of uniform phase. 

It follows from calculationJSI that for some model parameter values, 
the uniform phase becomes unstable with respect to fluctuations with q = 
(tt, tt) (the chess-board phase). 

To consider the thermodynamics of the chess-board phase analytically, 
we take into account the modulation of the pseudospin and electron dis- 
tribution, introducing two kinds of sites: (Sf) and ni correspond to one 
sublattice and (S"!) and n2 to the other one. In this case the single-electron 
Green's function for the I sublattice is equal to: 

Gi{k,c.„)^ (37) 

where gi{uJn) {I = 1, 2) is the nonperturbated Green's function for sublattice 
I. The single-electron spectrum is determined from the equation 
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The roots ei(ik) > £2(ik) > ^sitk) > £4(^k) of the equation ([55)1 form four 
subbands. The widths of subbands depend on the mean values of pseu- 
dospins. 

The branches ei(ik), £2{tk), on the one side , and e3(ik), £4{tw) on the 
other side, form two pairs of bands which are always separated by a gap. 
The equation for pseudospin mean values ([3T|l can be now written in the 
form: 

1 f « 1 _|_p-/3e ■ 

^Si;) ^ - tanh <! tL(h+nL-fyi ^ + In ■ 



where expressions for the effective self-consistent fields are 

2 * 



1 = 1,2 



(39) 



(40) 



4 = - 



e,(ik)+5(^f') 



(ei(ik)-ei(ik))(ei(<k)-ep(ik))(ei(ik)-£m(ik)) 
Expression for the electron mean number follows from ([33l 

4 



(ni+n2) = — ^ ^ r^[e»(tk)-/ 

k 4=1 

-2[((P+) + (P+))n(£) + ((Pf ) + (p-))n(£) 



(41) 



and the grand canonical potential ([27|l can be written for the two-sublattice 
case in the following analytic form: 



Nf3 



n cosh [f (£.(ik)-M)] 



E 

Z=l,2 

— In cosh I 
P I 2 



(cosh |£)2(cosh |e)2 

(/i+a^-a'i) +ln 
l+e-'^^ 



E(^')("2-«'i) 



-- Incosh |-(/i+4-a'i) + In ^^-^^ 



In 



(42) 



As previously in the investigation of equilibrium conditions, we consider 
two different thermodynamic regimes. 
a) The /x = const regime. 

The equilibrium is defined by the minimum condition of (|42p . Numeri- 
cal analysis of solutions of equations p9 |l - ([4T|) which satisfy this criterion. 
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was performed in Ref . . The examples of the calculated field dependen- 
cies of {SI—S2) (the order parameter for the chess-board phase) and grand 
canonical potential are presented in Fig. [T^] for low temperatures (the case 
g S> is considered) . It is seen from comparison of the $ values for the 
uniform and the chess-board phases that the modulated phase is thermody- 
namically stable at intermediate values of the h field in the region between 
points a and h. These points correspond to the first and second order phase 




0.25 0.30 h 0.25 0.30 h 

Fig. 12. Field dependence of the grand canonical potential and order parameter (T = 
0.005, = —0.36, W = 0.2, g = 1). Dotted and solid lines correspond to the uniform 
and the chess-board phases, respectively. See the text for a whole description. 

transitions, respectively. In the first case, the jump-like change of order pa- 
rameter is accompanied by characteristically similar changes of the subband 
widths and, as a result, by changes of the electron concentration. 

The resulting phase diagram [h — /z) at low temperatures is shown in 
Fig. [131 The chess-board phase exists as an intermediate phase between the 
uniform phases with different (S"^) and n values. Transitions between the 
uniform and modulated phases are of the first or second order and can be 
realized in the case when /i is placed in £2 or £3 subbands or between them. 
Transitions between different uniform phases (described in the previous 
section), which are of the first order, take place when the chemical potential 
is placed within the £1, £4 and partially within £2, £3 subbands. 

The shape of the (T—h) phase diagram strongly depends on the fi value. 
In the case when fj, is placed in the £3 subband, such a diagram is shown 
in Fig. I14h . With the temperature increase, the first order phase transition 
between the uniform and the chess-board phases transforms into the first 
order phase transition between uniform phases and, finally, disappears in 
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-g/2-W -g/2 -g/2+W g/2-W g/2 g/2+W 

Fig. 13. {h — fj.) phase diagram (T = 0.005, W = 0.2, g = 1). I: uniform phase, II: chess- 
board phase. Dashed hnes: first order phase transitions between the uniform phases with 
different pseudospin mean values. Dotted lines: first order phase transitions between the 
uniform and the chess-board phases. Solid lines: second order phase transitions. 

the critical point 9c. The diagram shows the possibihty of the first order 
phase transitions between uniform phases and either first or second order 
ones between the uniform and the chess-board phases at the change of 
temperature. 

In the case when chemical potential lies between the ei, £2 and £3, £4 
subbands, the transitions between phases I (uniform) and II (chess-board) 
are of the second order; the corresponding (T — h) diagram is shown in 
Fig. [Mb. 

b) The n — const regime. 

In this regime, the equilibrium is defined by a minimum of the free energy 
F = ^ + ^N. This condition forms a set of equations ([39| and (|4T|) for the 
pseudospin mean values and chemical potential. The obtained dependencies 
of F and fi on the electron concentration are presented in Fig. [T5l 

One can see the regions with d/i/dn < where the phase separation into 
regions with different phases (the uniform and the chess-board phases in 
this case) and with different electron concentrations and pseudospin mean 
values take place. 

Based on the obtained results, the phase diagram (n—h) was constructed 
(Fig. [16]). The phase separation into regions with the uniform and the chess- 



2, 2008 3:3 



WSPC/Trim Size: 9in x 6in for Review Volume 



stasyuk'new 



24 /. Stasyuk 




Fig. 14. (a): Phase diagram (T - h) (p = -0.36, W = 0.2, g = 1). L uniform piiase, 11: 
chess-board phase. Dashed lines: first order jjiiase transitions between different uniform 
phases (bistability) . Dotted line: first order phase transitions between the uniform and 
the chess-board phases. SoUd lines: second order phase transitions, (b): Phase diagram 
(T -h)iiJ. = -0.28, tk=o = 0.2, g = 1). 




Fig. 15. Dependence of the chemical potential fi on the electron concentration n and 
deviation of the free energy from linear dependence (T = 0.005, h = 0.28, W = 0.2, 
g = 1). Dashed line: phase separation area. Dotted line: uniform phase. Solid line: chess- 
board phase. 

board phases takes place when the chemical potential is placed within the 
subbands £2, £3. That agrees with the results obtained in the fi = const 
case when within this area we had the first order phase transition between 
the corresponding phases. 

The phase separated and the chess-board phase regions narrow with 
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n 



3.0 



1.5 



0.0 

0.0 0.2 1.8 2.0 

Fig. 16. Phase diagram (n — h) (T = 0.005, W = 0.2, g = 1). I: uniform phase, II: 
chess-board phase, PS: phase separation area. 

the temperature increase, but thick solid hues in Fig. [12] approach faster 
one another and, for high enough temperatures, we have only the phase 
separation into the regions with uniform phases. 

The fact of the existence of doubly- modulated (chess-board) phase in the 
PEM as intermediate one between segregated phases (at low temperatures) 
is similar to the results obtained in DMFT for the Falicov-Kimball model by 
Freericks and Lemansky.t^ Investigation of temperature Tc of instability of 
the high-temperature phase of the FK model as a function of the ordering 
wave vector (determined by the divergence of the relevant susceptibility) 
showed that near half-filling {pe = 0.5) the chess-board phase is stable, 
while the segregated phase exists at the occupations near electron band 
edges. The region of mo dulated phase narrows at the increase of g and 
disappears at o—^ oo.^^In this limit the thermodynamically stable states 
were analyzed^ based on the calculations of the grand canonical potential 
and the phase separation diagrams (where besides the spinodal lines the 
lines of the first order transition temperatures were shown) were built; the 
cases of different electron concentrations pe were considered. 

Nevertheless, one cannot claim the one-to-one correspondence. The cal- 
culations of Freericks and Lemansky were performed in another thermody- 
namic regime: the concentration of heavy particles pi was fixed (that would 
correspond to the given value of (5*^) for PEM). The closest to our case 
are the results obtained by Brandt and Mielsch.I^ They constructed the 
(T — Pi) diagrams at the finite g values for the FK model. The obtained 
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sequence of phases (segregated-modulated-segregated) when pi increases 
(from Pi = up to pi = 1) is the same as for PEM at the increase of the 
field h. 



3.2. Weak Coupling Case; U — 

Now let us consider the case of weak coupling {g < W). We can base on 
the perturbation theory approach taking the mean field Hamiltonian as the 
zero-order one. We use the approximation 

gn^S: ^ gn, (5f ) + g{n,)S^ - g{n,) (5f ) (43) 

in this case, based on the arguments that for the simplified PEM with U ^ 
the interaction constant g plays the role which is similar to that of U in the 
Hubbard model, and at the decrease of g below the critical value {g ~ W) 
the system should pass to the mean field regime of the Hartree-Fock type. 

Having in mind the possibility that the system can be in the uniform 
or modulated state, we shall consider, as above, two different cases. The 
first one corresponds to the homogeneous pseudospin ordering and spa- 
tially uniform mean distribution of electrons. The second one is the case of 
modulation with doubling of the lattice period (chess-board phase). 



3.2.1. Uniform Phase 

The Hamiltonian of PEM in the mean-field approximation (MFA) reads 
H = H,i + Hs + U, (44) 

Hel ^^{gV- Ai)"-i,<T + ^i^J^taO'O.'^ ' 

i.(7 'i'^3,(y 

H,=Y.^{gn-h)St + nSf\, (45) 

■i 

U = — g nrj = —Ngnrj. 

i 

The parameters (ui) — n and (5f ) — rj are determined from the set of 
equations 



1 ^(e^(..+..-M) + l)-i ^ i- ^ figrj + t^), (46) 



N ' N 

kcr ko 



rj^^^tanhi^): X = ^{gn - + Q^. 
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The grand canonical potential in the MFA is given by the expression 

^ = - 1 ^ ln(l + c"^^ ) - T ln(2 cosh ^) ~ gnr^ (47) 

k.cr 



where A — y/ {gn — h)^ + il^. 

Similarly to the strong coupling case, we can distinguish the regimes 
of the constant electron chemical potential /i = const (where the stable 
states can be found from the minimum $ condition) and the given electron 
concentration n = const (when one should find the minimum of the free 
energy F = /in + 

a) Thermodynamics in the n = const regime. 

We consider at first the simplest case fj, = 0, h — g. Here the solution rj — 
0, n = 1 of the set of equations exists at any temperature and describes 
a disordered phase. Furthermore, at low temperatures there appears a non- 
zero solution 77 ^ 0, n ^ 1. A critical temperature Tc is determined from 
the equation 

l + ^tanh^no = 0, (48) 
which in the case J7 ^ reduces to the form 



Here 

Ho = 



1 + f^g'Uo = 0, (49) 



r, rW 

^/'(tk)=2/ dtp{t)f{t). (50) 
J-w 



k 

2 

In the low temperature limit we have Tc = ^p(O) (when the DOS at 
the Fermi level p(m)|^_o is finite). In the case 7^ there exists such a 

critical value flcr — 5^p(0) {^cr = 2w rectangular DOS), above 

which (at H. > ficr) the phase transition to ordered phase disappears. This 
is equivalent to the existence of a critical value of at given the phase 
transition is possible when g > gcr — \J '^yj- *^^se of the DOS with 

logarithmic singularity, the critical temperature exists at any values of the 
tunneling splitting parameter and at ^ 3> we have an asymptotic 
expresssion: 

r,«2el^exp(^^^). (51) 

The physical nature of the phase transition considered here aX h — g^ ^ — Q 
(the fixed jjt regime) is as follows: the appearance of an ordered phase is 
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connected with its stabilization due to the shift of the electron band down 
to the low energy values under the effect of the internal field; this ensures 
the corresponding gain in the electron energy (it should be mentioned that 
the electron band spectrum in this case remains unsplit in the uniform 
phase; only the shift of the band as a whole can take place). 

This mechanism remains the main reason of the phase transition when 
the initial electron band is not half-filled. In this case (when 7^ 0) we 
performed the investigation using the numerical calculations when the set 
of equations (|46p is solved and using the expression (|47p for the grand 
canonical potential The selection of solutions was carried out using the 
condition of the a bsol ute minimum of <&. 

As can be seenp^l below the system undergoes the first order phase 
transition with jumps of the mean values of the electron concentration 
and pseudospin at the change of the field h] the phase transition point is 
determined using the Maxwell rule. The similar transition takes place at the 
change of the chemical potential at fixed h. The presence of the tunneling- 
like splitting decreases the temperature of the phase transition at the fixed 
values of /i and h. 

The regions of coexistence of phases with different values of the electron 
concentration and pseudospin are shown in the plane {fi,h) at f2 = and 
n^Oin Figs, [m 

Phase transition lines in the plane (T, h) at different values of n are 
shown in Fig. [181 Such a line is vertical for the case ^ = only; for the 
case /i 7^ 0, the line is bent. This makes the first order phase transition 
possible at the change of temperature (with the jumps of the parameters 
r], n). The slopes of the phase equilibrium curves are opposite for /x > and 
/i < 0. The lines of the critical points arc shown for the cases = 0, $7 7^ 0; 
the calculations are carried out using rectangular density of states. Similar 
phase diagrams are obtained using direct momentum summation for square 
lattice (Fig. [T9|) . As in the case of rectangular DOS, the maximum Tc is 
achieved at /j, = but the critical temperature line has a more pronounced 
peak. 

b) Thermodynamics in the n = const regime. 

The values ni,n2 and ?yi,?72 between which the jumps of the electron 
concentration and pseudospin mean value take place at the phase tran- 
sitions in the fi = const case, correspond to the phases which coexist in 
the phase transition points. In the regime n = const, there is a phase 
separation on the phases with the above mentioned values of n and fi. 
For example, at the parameter value H, = 0, h = 0.7, T — 0.008, the 
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-W-ig/2 -W+g/2 
(0 \ -,/2 



Fig. 17. (a): The ground state diagram (T = 0, f2 = 0). Regions with different n,rj 
values are separated by dashed lines and a solid line (the phase transition line); (1) 
n = 0,»? = -l/2; (2) n = 0,7? = 1/2; (3) r, = 1/2, n = 1 + - 5^; (4) r, = -1/2, n = 
1 + + (5) n = 2,r) = -1/2; (6) n = 2,r; = 1/2. (b): The {h - ^i) phase 

diagram (T = 0.004, f2 = 0.12). Regions with different n,r] values are separated by 
dashed lines and a solid line (the phase transition line): (1) n 0, ?? ~ h — 1!-^ = ; (2) 



nKi 2,r} ; 



h~2g 



2 VCi-ajj^+n^ • 

0.07 



0.00 



-0.8 


1" 


=0 










n=o.i 



0.10 



0.50 
h 



At=0.8 



1.10 



Fig. 18. The critical temperature lines (dashed lines) and the phase transition lines 
(solid lines) for Q, = and H = 0.1 (the case of the rectangular DOS). 



system is unstable with respect to the phase separation in the region 
m = 1.149 < n < ?i2 = 1.645.1221 
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Fig. 19. The critical temperature lines on the a) {h, T) and b) {fi, T) planes for H = 
(direct momentum summation is used). 

Phase separation regions are shown in Fig. [20] at different temperatures 
(the calculations were carried out for square lattice using direct momentum 
summation when solving the set of equations (|46p ). At the increase of tem- 
perature the separation region narrows and at T > Tc it disappears. The 
presence of the tunneling-like splitting leads to the decrease of the area of 
phase separation region and to the lowering of Tc- 

It can be noted that at weak coupling, when only uniform states are 
considered, we have one phase separation region in the (/i, n) plane while 
in the strong coupling case (17 ^ W) there exist two such regions at the 
distance of order g along h axis; the modulated phase is placed between 




Fig. 20. (n — /i)-phase diagram in the cases a) SI = 0; b) SI = 0.1. Phase separation 
regions for different temperatures are shown, a) 1: T = 0.008, 2: T = 0.08, 3: T = 0.0972; 
b) 1: T = 0.008, 2: T = 0.075, 3: T = 0.0903. 
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thenJSH (see also Sec. [5X^ . 



3.2.2. Phase with Double Modulation 

Let us consider now the thermodynamics of the simphfied PEM in the case 
of doubly modulated phase. The possibility for such a phase was shown at 
a strong coupling {g ^ W); the necessary condition was the location of the 
chemical potential between the split electron subbands. Now (at g < W) the 
band is unsplit and in such a situation there must be another mechanism 
of stabilizing the lattice modulation. 

At the double modulation the crystal can be divided into two sublattices 
(a = 1, 2), and the parameters rja = (S^^^), Ua = '^{niaa) can be introduced 

{i is an unit cell index). Similarly to the homogeneous phase, the mean field 
approximation is used. The modulation leads to the splitting in the electron 
spectrum due to difference between the internal field acting in sublattices 

Aka=5 — ^ + (-1)y(5 — ^ ) +*k (52) 

(the similar effect takes place in the FK modeP^. The initial band is 
divided into two subbands separated by the gap A = g\rji — 772] (see Fig. 

EB). 



band edges 

w [(i/2g(Ti,-n,)f+w'i" 



i/2g|ii,-ij 

Q 

1/2g|Ti,-nJ 



-W ■ 



-[(1/2g(ri,-il2)f+W^]"' 



T>T T<T 



Fig. 21. The band edges for the cases T > T* (homogeneous phase) and T < T* (the 
double modulation case); fi = 0, h = g. 



Contributions from both sublattices are present in the equation for the 
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electron concentration in sublattices 



N 2 

ko- 



(53) 



where 



cos 20 = 



y 2 



can be obtained diagonalizing the mean-field two-sublattice Hamiltonian of 
the model. ^"^ Another equation which appears as a result of the averaging 
of the operator S'f^ has the form 

,. = ^tanh(^). (54) 

Here 



K = Vianc - hf + . (55) 

In the mean field approximation, the grand canonical potential for the 
double modulation case has the form 

^ = -^2]ln((l+e - )(l + e - )) 

k,(T 

— Tln(4cosh cosh ^^) ~ dinirji + 77.2772)- (56) 

The solution of the set of equations for the Uq, and 77^ parameters 
and the inv esti gations of thermodynamically stable states were carried out 
numerically!^ It was established that phase transitions from the uniform 
to the low temperature modulated phase are of the second or first order. 
This is illustrated in Fig. [511 where the phase transition lines at /i = are 
shown for f2 = and 7^ (here a direct momentum summation is used 
in calculations); the tricritical points are also present here. 

The difference Sn = ni — n2 (as well as the difference 5rj = rji — 772) 
can play a role of the order parameter for the modulated phase. Coming 
from the equations for 5n and St] we obtain the following condition of the 
appearance of nonzero solutions 

kr 



(57) 
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Fig. 22. The phase transition hnes (sohd and dashed hnes are the hnes of the second 
and of the first order phase transitions, respectively) from the uniform phase to the phase 
virith double modulation (1 : f2 = 0; 2 : C = 0.2). 



Proceeding from this equation, we can find a critical temperature T* as the 
maximum temperature (among the set of temperatures which are obtained 
for different h values) which fulfills this equation at a fixed value of the 
chemical potential. This temperature is the point of the second order phase 
transition to modulated phase at the corresponding value of the field h. 

In the symmetric case, when fi — 0, h — g and in the high-temperature 
phase n = 1, 77 = 0, the equation (|57p reduces to the form 

1— S*-i^?^E^/(*'^)- (58) 



The critical temperatures T*, obtained in the cases of rectangular DOS and 
DOS with logarithmic singularity, are, respectively, 

eW 



■ exp 



9^ 



and 



T* « 2eW^exp 



9 



(59) 



(60) 



In both cases, they are higher than the corresponding temperatures Tc for 
transitions between uniform phases and remain finite at high n values. 

It can be seen that T* > Tc also at /x ^ 0, but for the fj, values, which 
are less than the certain (|/i| < /xq) value. The typical dependencies of T* 
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and Tc on /i are shown in Fig. [221 in the cases = and ^ 0. One can 
conclude that in the case of the electron band occupation close to the half- 
filling, the transition to modulated phase should be realized. The transition 
between two different uniform phases is possible only when /i is placed near 
the band edges. 




(a) (b) 

Fig. 23. The dependence of the critical temperature on the chemical potential; a) f2 = 
b) O = 0.1. Solid line refers to the case of the phase with double modulation, dashed 
line refers to the transition into the homogeneous phase. 

An existence of the first order phase transitions between uniform and 
doubly modulated phases, shows the possibility of a separation into these 
two phases. This takes place at certain values of the electron concentration. 
The corresponding {n — h) diagrams are shown in Fig. 1241 The borders of 
the separation regions were obtained from the convexity condition of the 
free energy defined as F/N = ^/N + n^. As temperature increases, the 
separation area narrows, but in the middle of it there appears a region of 
the chess-board phase existence. This is an additional feature which supple- 
ments the picture of separation shown in Fig. 1201 At high enough temper- 
atures, only the second order phase transition into the doubly modulated 
state remains and the phase separation region disappears. 

3.2.3. Pair Correlation Function and Susceptibilities 

The analysis of thermodynamically stable equilibrium states of the PEM in 
the case of weak coupling can be supplemented by an investigation of tem- 
perature and wave vector dependencies of the pseudospin, electron density 
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0,3 0,4 0,5 0,6 0,7 ' 0,3 0,4 0,5 0,6 0,7 

h h 



Fig. 24. (ra— h) phase diagram, Q = 0. Phase separation regions are shown for different 
temperatures: a) T = 0.008, b) T = 0.08. Dashed lines denote the borders of the uniform 
phase, thick solid lines denote the borders of the phase with doubly modulated lattice 
period. 



and mixed pair correlation function. The corresponding Green's functions 
were calculated in Ref. 33J within the GRPA scheme. The cases of isother- 
mal response were considered, where the isothermal susceptibility 

XT(q,WrO= / (T,M(0)M(r))qe^'^"Mr - /3(Af)2,5K) (61) 
Jo 

is expressed in terms of the Matsubara Green's function and the so-called 
"isolated" response, which is described by means of the two-time Zubarev 
Green's function 

X/(q,o.)^ ((M|M))q^. (62) 

The dipole moment Mi of the unit cell was taken in the form: Mi = 
deUi + dgS^; here the electron contribution due to nonhomeopolarity of 
occupancy of the electron orbitals was taken into account besides the pseu- 
dospin contribution. 

Such an expression for dipole moment comes from the form of trans- 
verse component of polarization in the case of the YBaCuO structure. The 
electron component corresponds to the charge transfer in a perpendicular 
direction from/to the CU2O2 layers (with the participation of the Cu-0 
chains), while the pseudospin component is connected with the redistribu- 
tion of the ionic charges when O4 ion moves from one equilibrium position 
to another.lZHES] 

Similarly to the above considered strong coupling case, the simple se- 
quences of loop diagrams are taken into account in the diagrammatic rep- 
resentation for Matsubara's correlators. The connections between loops 
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are accomplished by semi-invariants or by the boson (pseudospin) Green's 
functions. The contribution that corresponds to the separate link is (see 
Ref. 

Sq(cc;) = \ sin^ ^[Kl{ujn) + if^(-t^„)](a")o - Afq(t^„). (63) 



Here 



^qK) = -■ ^^qK) = Pb' cos" d5{LUn), 

sm^^n/X; (a^)o = 6=itanh^; = (64) 

The pseudospin Green's function K^{ujn) is constructed of operators of 
the transverse pseudospin components acting in the rotated reference sys- 
tem {K° - (Tct+ct-); CTf = S*! cosi9 - 5f sini?; af = Sfcos^ + S'f sini?) 
and has a pole at A (the pseudospin reversal energy). The semi-invariant 
Mq{ujn), describing the correlation of longitudinal pseudospin components 
(~ T{a^a^Y), is proportional to 6{uj). There is no such contribution 
in Sq(ti') when the Green's functions are calculated using the equation 
of motion and the decoupling procedure for Zubarev's functions. Thus, 
the isothermal (xt) and isolated (x/) susceptibilities do not coincide for 
PEM.I33] should be mentioned that a similar result is also obtained in the 
cases U and 



In the GRPA scheme, the summation of loop sequences leads to the 
expression for pseudospin correlator 



l-<?2l]q(,^)nq(cc>)' 



Eq(w) 
^q( 

2 ^ n(tk) - "^(ik-q) 



where 



nq(a;) = ^ > ' ' (66) 

is the standard electron loop contribution. 

The condition {TS^ S^)q^uj=o oo (which corresponds to divergence of 
isothermal susceptibility xt) indicates an instability with respect to tran- 
sition into modulated (at q =/= 0) or another uniform (at q = 0) phase. 
The thermodynamic parameter values, at which xt ~* oo, determine the 
spinodal points. 

The equation 

A^ + 5^ sin^ 6'A(a^)onq + X^g^pb' cos^ 9Uq = (67) 
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was solved together with equation (|46p for the mean values 77 and n writ- 
ten in MFA. The function Ilq was calculated numerically by the direct 
momentum summation for square lattice. 

At fixed values of chemical potential, the critical point can be defined as 
an upper point of spinodal (on the (T, h) plane) with the highest tempera- 
ture depending on the wave vector q value. Fig. [25l shows the dependencies 
of the critical temperature and the corresponding wave vector on the chemi- 
cal potential in the case fi = (only positive values of the chemical potential 
are shown; at h = g, the picture is symmetrical with respect to the point 
/i = which coincides with the centre of the energy band). We can see that 
the case q — (tTjTt) is realized when |^| < 0.25 at chosen parameter values, 
which means that the system can pass into the phase with doubly mod- 
ulated lattice period. The case q = (transition into the uniform phase) 
is realized when 0.85 < < 1.250 The system undergoes transition to 
the incommensurate phase at intermediate values of the chemical poten- 
tial. The presence of tunneling splitting narrows the interval of values of fx 
at which the above mentioned transitions take place; at high enough values 
of Q, the transition into the chess-board phase occurs only. 




0.00 0.40 o.so J.20 0.00 0.40 o.so J.20 



Fig. 25. The dependence of the modulation wave vector q=(q,q) and the temperature of 
absolute instability of high— temperature phase on the chemical potential, Q = 0, g = 0.5. 

The given results a,t — generally correspond to the picture of 
phase tr ansitions in the FK model obtained in DMFT in the case of weak 
couphng.EZESIMlit was shown for FK model that at small values of n the 



^1.25 = W + g/2; this value corresponds to the upper edge of the band when (S^) = ^. 
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phase separation can be realized, while near half-filling of the band the 
chess-board phase is preferable and, finally, at intermediate values of n the 
appearance of phase with an incommensurate modulation is possible. The 
transitions into one or another phase were observed by the divergences of 
corresponding susceptibilities. It should be mentioned that such a proce- 
dure did not enable the authors to reveal the thermodynamically stable 
states in the regions where the instabilities of the both types (at q = and 
q = (7r/a,7r/a, ...)) are superimposed; this problem can be solved based on 
the analysis of the behaviour of the grand canonical potential. 

The GRPA scheme used here to investigate the PEM is advantageous 
in interpreting the dielectric susceptibility divergences due to the explicit 
dependence of the XT(q, 0) function on the wave vector. In the DMFT 
approach at c? — *■ cxd such a dependence enters only through the function 

d 

X{q) = 3 X] cosgj that leads to some difficulties in considering the incom- 

mensurate ordering. Besides, in Refs. [571 1631 166j there was used a regime 
of a fixed concentration of localized particles; in the PEM this corresponds 
to the regime {S^) = const. In this case, the authors came to a conclusion 
that the transition to the chess-board phase at 17 < VF is always continuous 
and spinodals are the lines of phase transitions.^^ In such a situation the 
separation into the uniform and modulated phases would be impossible. 
Contrary to that, by analysing the behaviour of the grand canonical poten- 
tial we showed that the transition to the chess-board phase can be both of 
the second order (the spinodals are the phase transition lines) and of the 
first order (spinodals do not coincide with phase transition lines). Due to 
the first order phase transition to the chess-board phase, the possibility of 
the phase separation into the uniform and the chess-board phases in the 
case of weak coupling was demonstrated. Such a possibility, as was shown 
in Ref. [M] , see also Sec. 13.1.21 exists in the case of large values of coupling 
constant. It should be mentioned that in Ref. [BT| only the possibility of 
phase separation into different uniform phases was investigated. 

Though the phase transitions in the PEM at weak coupling are similar 
to the transitions revealed in this model in the case of strong interaction, 
g ^ VF,^*^ the physical mechanisms of transitions are to a greater extent 
distinct. In the case of strong coupling, the electron spectrum is always split 
due to the one-site interaction. The mechanism which ensures the advantage 
of the transition is connected with the different character of the electron 
spectrum reconstruction in the subbands and with the corresponding redis- 
tribution of the electron density of states. At weak coupling, a new phase. 
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which appears at the transition between uniform phases, is stabihzed due 
to the shift of the electron band as a whole. The phase with a double mod- 
ulation appears due to energy gain at the splitting of the initial band at the 
Fermi level (the effect is similar to the Peierls instability at the interaction 
with phonons). Besides that, the dependencies of the critical temperatures 
on the coupling constant g are different in the both cases: (or T*) is 
proportional to a.t g ^ W, while a,t g ^ W the critical temperatures 
decrease (~ i) when g increases (such a type of behaviour of T* for the 
FK model was obtained in Refs. [5T1 163| ). 

3.3. Superconductivity in the PEM 

Due to the presence of intrinsic dynamics, the PEM at il 7^ also possesses 
an instability with respect to transition into superconducting (SC) state. 
Under certain conditions the transition to SC state will compete with the 
transition to modulated phase (CDW). Among others, close attention to 
this problem for an electron system interacting with anharmonic structure 
units was paid in Ref. [67j . It was shown that, when there is no electron 
correlation, the transition temperature to CDW, Tp, is higher than that 
to SC, Tc- The authors, however, considered only the case where the lo- 
cal potential is symmetric and the electron filling is close to a half. The 
question concerning the appearance of superconductivity for a wider range 
of the model parameters and electron concentrations has not been exam- 
ined. It was just the point of investigation performed in Ref. |35 within 
the GRPA. To examine the possibility of the appearance of SC phase we 
calculated the static susceptibility x^'~^ i^i the superconducting channel. 
Taking into account the diagrams which correspond to ladder approxima- 
tion (with the parallel directions of lines of the fermion Green's functions), 
we obtained the Bethe-Salpeter equation for the superconducting vertex 
part r^,,^,(ki,k2|2l 

r = + T ^ r°x°r, (68) 

where xlS^i) = j^Gl^{^i)G\^[~LOi)- = g^{TS^S^)u,-k^.u,-.., is 
given by the expression (j65p . The susceptibility 

^^"^ = E / (T.axk(r)«T-k(r)4_qa+ )e-"-dr (69) 
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is connected with the vertex part F in the following way 

^X^^ = Ex-W+^ E X°,(ki)r.„.,(ki,k2)x°,(k2). (70) 

aj,k ki ,k2 ,0^2 

Both approximations, the ladder one for the Fi^j.^jj (ki, k2) vertex con- 
struction and the chain (GRPA) one for the construction of ^2(ki,k2), 
were employed by analogy with what was done in the t — J model for 
the SC description!^ The mentioned approximations correspond to the 
known Migdal-Eliashberg (ME) one, which is usually used in considering 
the electron-phonon systems, particularly in the Holstein model. Consid- 
ering the phon on frequencies to be small in comparison with the trans- 
fer integral!^ one gets results in qualitative agreement with the quantum 
Monte Carlo simulations. Since the PEM is similar to the Holstein model, 
and may be considered as a double-level approximation of the latter one, the 
ME approximation is expected to be satisfactory in the limits of non-half 
filling, low temperatures, and small 

A way to find the SC transition temperature is to determine a temper- 
ature, at which the susceptibility in the superconducting channel diverges. 
It corresponds to the condition when the scattering matrix 

T^i^2(ki,k2) = rx;i,(ki)r;i^„^(ki,k2) (71) 

has an eigenvalue which is equal to unity.^^^^^ In calculations performed in 
Ref. 134] ■ the approximation « —g^Ti, sim ilar to the non-renormalized 
ME approximation in the Holstein modelj^ was used. This can be done 
in the high-temperature phase, when the system has not still passed to the 
CDW state. In such a case the unperturbed vertex part does not depend 
on the wave vector and thus the unit eigenvalue of the matrix 

TL,c.2 =7^Ex°2(ki)r°,.2 (72) 

ki 

should be found at first. Of all the temperatures within the (T, h) plane 
which satisfy this condition, the highest one is chosen as the critical SC 
transition temperature. 

Numerical calculations performed in Ref. [34] show that the SC tran- 
sition in the simplified PEM is possible (in the weak coupling case) at 
the electron concentrations away from half-filling (when the chemical po- 
tential is placed near the electron band edges) and outside the region of 
the above described transitions with the modulation of the electron and 
pseudospin density (see Fig. [26| . Such a picture is similar to that obtained 
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within the DMFT for the Holstein modef^ (but, in the latter case, the 
incommensurate phase does not appear at the intermediate values of ^). 
This is also consistent with the results obtained by quantum Monte-Carlo 
simulations P^IZll These papers just established that SC could appear at low 
temperatures; the transition to incommensurate CDW was not studied. 




Fig. 26. The temperature of the absolute instability of the high-temperature phase with 
respect to the transitions to CDW (solid line) and SC (bold dashed line) as functions of 
the chemical potential, for Q = 0.2 and g = 0.5. Dotted line shows the temperature of 
the transition to the phase with the modulation wave vector q = (tt, tt). 



The estimates performed in Ref. show that at « 0.5 eV and g « 
0.25 eV, the maximum value of the SC transition temperature is T^'^^ « 
10 — 40 K. This agrees with the conclusions obtained for PEM in Ref. [57] . 
where calculations were performed by analogy with the scheme used for the 
Eliashberg equations in the limit of a weak electron-phonon interaction, 
provided that the renormalization of the pseudospin excitation energy is 
neglected. According to the estimations carried out in Ref. [72], if there 
were no CDW, the SC transition temperature would be T^^ « 40 K at the 
band filling close to a half. 

The question of what occurs in the systems as the tunneling splitting 
frequency grows further is also of great interest. It was shown in Ref. |32| 
that, at /Lt = 0, the critical temperature of transition to the phase with 
double lattice period modulation, T^^"^ , decreases with the increase in 
n according to the exponential law: Tf^^ - exp /'- ^vTw ^ see §0^. It 
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follows from the analysis of the behaviour of the T^^^i^^ matrix elements^ 
that the SC transition temperature, T^*^ ^ changes with in a similar way. 
In this case, for ^JL = 0, T^'^ {Q oo) ^ T^dw (f] ^ qo). As is seen 
from Fig. I26[ with the increase in \fi\, the CDW transition temperature 
falls more rapidly than the SC one. Thus, for the non-zero values of the 
chemical potential, provided that is sufhciently large, the SC transition 
temperature is expected to be higher than the CDW one, and there will 
only be the transition to SC. However, as was noted above, to make the 
correct analysis of the competition between these transitions, the renor- 
malized F*^ vertex should be used when solving Eq. (j7ip and determining 
Tsc- Moreover, when fl values are sufficiently large, the applicability of the 
approximation in deriving this equation turns out to be unjustified. 

The mechanism that leads in the PEM to SC, which we have considered 
here, as well as the traditional phonon one, does not result in high values 
of Tc, and apparently it does not explain the HTSC phenomenon. 

4. Thermodynamics of PEM at Finite U Values; the 
U — >■ oo Limit 

The analysis similar to the one given above was performed in Refs. |24| - 
I26l73tf76] for the PEM with U ^ 0. The presence of the electron-electron 
on-site interaction leads to some differences in the behaviour of the model 
with respect to the case of FK model (even if = 0). The consideration 
was based on expansions in terms of electron transfer tij, as in the strong 
coupling limit g 3> VF. 

The single-site Hamiltonian 



Hi = Uni^riii - fi{ni^ + "ij) + g("-iT + nil)Si - hS- - ilS^ (73) 



that includes the U -term can be redu ced to the diagonal form using the 
rotation transformation \R) = a^^^^ \r)^^^^ where \R) — \ni-^ ,nii, S^) is 
the single-site basis of states 



|1> 
|2> 
|3) 
|4> 



10,0,1/2) 
11,1,1/2) 
10,1,1/2) 
|l,0,l/2) 



|1) 
|2) 
|3) 
|4) 



|0,0, 

ll,l, 

|0,1, 

|i,o. 



1/2) 
1/2) 
1/2) 
1/2), 



(74) 



and 



rirg — h 



(75) 



cos (fir 
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In terms of Hubbard operators X'^'^ — \r){s\, acting on the new basis, the 
transformed Hamihonian of the model is 

ir ij a 

with 
and 

rnn mn 

Here ni = 0, n2 — 2, = = 1 [uf = rir); expressions for are given 
in Ref. [26]. 

The detailed consideration of single-electron spectrum of model ([T]) at 
U ^ was performed in Refs. [23l [74] in the Hubbard-I approximation. 
The interaction with the anharmonic (pseudospin) mode splits the energy 
levels of ordinary Hubbard model 0, Eq and 2i?o + U into sublevels A^^f . 
As a result, each Hubbard single-electron band splits into four subbands. 
In the independent subband approximation 

£,,(q) = A, - A, + tq {A-,,)^{X'- + X''), (78) 

where (rs) = (41), (41), (41), (41), (23), (23), (23), (23) for a =t and 
substitution 4^3 should be done for a =|. The widths and statistical 
weights of the subbands are determined by parameters A^^. The four lower 
and four higher bands correspond to the hole and electron pair motion, 
respectively. Their positions depend on the value of the asymmetry field h 
(see Fig. [27] where the band edges as functions of h are shown at the fixed 
value of /x) . 

Thermodynamics and dielectric properties of PEM with U ^ were 
studied in Refs. [26l [75] . Being interested in the problem of lattice insta- 
bilities in high-Tc superconductors of YBaCuO - type we calculated the 
transverse dielectric susceptibility (corresponding to the Szz component 
for YBaCuO - structure), using, as above (see Sec. I3.2.3p . the expression 
Pi = dsSf + deUi for local dipole moment. We can respectively separate 
ion, electron and mixed components in the total susceptibility 

X±(q,a;„) = 4x'^'^(q, ^^n) + rfeX""(q, t^n) 

+ dsd, (x^"(q,c^„) +x"^(q,'^n)) , (79) 
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Ers{±W),ll 




Fig. 27. Single-electron spectrum depending on h/g in the fi =const regime; U = 2.2, 
3 = 1, h = 0.7, Q = 0.3, W = 0.2, ^l = 0.5. 



where x — K in the regime /i = const (when we fix the value of the 
chemical potential and permit the charge redistribution between conducting 
sheets CU2O2 and other structural elements, having in mind an application 
of the model to the YBaCuO-type crystals) and 



lA' 



(80) 



in the regime n = const ( whe n we fix the electron concentration n in the 
conducting sheets Cu202).^Here K^^ (q, ^n) are Fourier-transforms of 
semi-invariant Matsubara's Green's functions 



(81) 



constructed of the operators , rii. 

In the case of zero hopping [tij — 0) the exact expressions for correlation 
functions x^'^i'^n), X^^i^n), X^^i'^n) can be easily obtained.^ 

For n=const regime the main contribution into susceptibility is pro- 
duced by the pseudospin subsystem. The susceptibility 



vc 



(82) 
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as a function of h generally possesses three peaks. The maxima at /i = 0, 
g and 2g correspond to the points of possible dielectric instabilities (that 
can appear due to the proximity of corresponding energy subbands). The 
intensities of peaks are redistributed with the change of n; this is illustrated 
in Fig. EHli for the case U/W > 1 (fi = 0.3^, < n < 1), when only two 
peaks are present; the first peak disappears at n ^ 1 and the second one 
disappears at n — > 0. The effect remains the same at Uj —^ 0. 




(a) (b) 

Fig. 28. (a): Dielectric susceptibility as a function of h/g at different electron concentra- 
tions; U = 2.2, g = 1, h = 0.7, n = 0.3, W = 0.2, ds = 0.6. (b): Pseudospin component 
of susceptibility Xa9 vs h/g for r2 = {uj„ = 0, n = 0.9, T/g = 0.05). 

In the /i =const regime there are three regions of h values (h < 0, 

< h < g and h > g) with different field and temperature behaviour of 
susceptibility (see Fig. [28b). This is a result of the difference in the direction 
and value of the total effective field h — Uig acting on the pseudospin at site 

1 when electron occupation on this site changes (n^ = or 1), which testifies 
to a strong correlation between electron and pseudospin subsystem. 

The behaviour of transverse dielectric susceptibility x"^ as a function 
of the field h and temperature becomes more complicated when we take 
into account the electron tran sfer. T he calculations of components of sus- 
ceptibility ([79|) were performecPSESl a,t tij ^ in the particular case of no 
tunneling splitting in anharmonic potential well (17 = 0). 

Green's functions (|8T|) . constructed of the operators 

4 4 

n,,=^n.(Xr + Xn, St ^ ^n^iXr ~ XD, (83) 

r— 1 r— 1 
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can be expressed in terms of functions 



{TXf''{r)X^J{r')am, 



(84) 



The perturbation theory with respect to electron hopping term ij,- and 
the corresponding diagrammatic technique for Hubbard operators^ were 
used to calculate these functions. The diagrammatic series were summed 
up in the GRPA, where the sequences of electron "loops" connected by 
semiinvariants (similarly to the case of the simplified PEM, see Sec. I3.1.1|) 
or by vertices with the three fermion lines (such vertices appear only at 
U ^ 0) are taken into account. 

The Fo urier t ransforms K^P'^\q,u)n) of fimctions (|5T|) can be presented 
in the fornP^ 



KPP (t^„,q) =/3|6pp,(q) 

+ pn(o,q)S(q)n'(o,q) 



'n(0,q)6(q) 



pp' 



6(q)n'(0,q) 



pp' 



pp' 



(5 (tj„) + n;'p, {uJn,^) 



(85) 



The first term in the right hand side of (l85|) is the "full" semi-invariant of 
the second order that satisfies the Dyson-type equation 



bllb 



(86) 



Here bpq = {XppX'^'^)qc is the second-order semi-invariant calculated in 
mean-field approximation and the full "loop" contributions H, H', 'H, H" 
are determined from the Bethe-Salpeter type equations 



H" 
H' 



n[,' 
n' 



n;,'H' 



'n = 'Ho 

H = no- 



'HoH", 
n'„ H' + Ho H", 
-'Ho 'H + Ho H", 
Ho 'H + H[, H . 



(87) 



Zero-order polarization loops Hq, Hq, 'Hq, Hq are constructed of the 
single-electron Green's functions which are calculated in the Hubbard- 
I approxirnation (corresponding to the summation of chain fragments of 
diagrams). 123261 Loops HJ,, 'Ho, HJ,' are determined only by intraband tran- 
sitions whereas loop Hq is determined by the interband transitions as well 

jj / ^ 1 >r^, , n+ [emr (k)] - n+ [e^p (k + q)] 



where £mr(k) are determined in ([78|l and ri+(A) is Fermi distribution. 
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Numerical calculations of the static dielectric susceptibility x± (Qi 0) 
performed in Refs. [25l [75] along the (0, 0) (tt, tt) line in the 2D Bril- 
loine zone, revealed that, similarly to the above considered case U — 0, the 
essential feature of the model is the presence of divergences on the tem- 
perature dependencies of functions Xn ^^'^ Xft (Fig- ES]) ■ They appear in 
a certain range of the model parameter values. For h > and h < g such 




0.01 0.1 1 



Fig. 29. Temperature dependence of x"" at n = 0.95 and q = ({/ — > oo, h/g = 1.05, 
W/g = 0.2 ). 

divergencies exist only at q = (F-point) and can be treated as the mani- 
festation of the dielectric type instabilities which appear in the pseudospin 
subsystem (i.e., the system of anharmonic oscillators) under the influence 
of the effective interactions. As in the case of a simplified model, it means 
that the system can transform to another uniform phase. The corresponding 
phase diagram T* vs n is shown in Fig. I30r a). where the lines limiting the 
stability region are plotted. For < h < g, besides the dielectric instability 
at F-point at n > 0, the instability at q = (■^, ^) (Af -point) with respect 
to the charge ordering (double modulation) occurs (Fig.lSpyb)). Such an in- 
stability is realized when the electron band (that corresponds to the lower 
Hubbard subband at [/ — s- cx)) is nearly fully occupied. 

The above mentioned divergences on the temperature dependencies and 
an increase of the dielectric susceptibility which takes place in the vicinity 
of h = and h = g values (for U — > oo, n < 1) are connected, first of 
all, with the polarization contributions Ho from the electron transitions 
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(a) (b) 

Fig. 30. Dependence of the instability temperature T* on the electron concentration 
n in the cases h > g and < h < g. (a) U oo, f2 = 0, n/g = 1.05 (dielectric 
instability, q = 0). (b) U — > oo, H = 0, n/g = 0.95 (o o o — dielectric instability, q = 0; 
* * *— instability with respect to CDW ordering, q = (it, it); dashed line shows the MFA 
results). 



between hole subbands (41) and (41) (or (41) and (41)). Just for h ^ and 
h ^ g these subbands come nearer to each other. 

It should be noted that the results presented here for dielectric response 
of PEM in the C/ — > cx) limit were used in Ref. [25] in describing the dielec- 
tric anomalies which have been observed in the high- Tc superconductors 
of the YBaCuO type at an early stage of investigations. In particular, the 
temperature dependence of Xn given in Fig. [29] was related to the one ob- 
tained experimentally for YBa2Cu307-^^ where the similar behaviour of 
dielectric permittivity Sc was found at the temperature near (and above) 
the point of the superconducting phase transition. 

An interesting feature of the model at large values of U is the possibility 
of change of the electron concentration n under the influence of h. Such an 
effect can take place in the /i=const regime; this becomes possible when the 
chemical potential leaves (or enters) the energy subbancP^ (Figs. [27] and 
IBlT a)). The change of n is accompanied by the corresponding change of the 
pseudospin mean value (Fig. [3lT b)). In the case of YBaCuO structure, this 
means that if fi moves inside the subband with the change of e.g. external 
electric field, then the electron concentration (average occupancy of states) 
in CU2O2 layers changes. This causes a redistribution of O4 ions in their 
equilibrium positions. All this may correspond to the so-called electric-field 
effect observed in HTSC compounds (see, for example, Ref. [77]). 
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(a) 



(b) 



Fig. 31. Mean electron concentration n (a) and pseudospin mean value (b) as a function 
o{ h/g in the ^=const regime. 



The obtained results show that in the U ^ oo hmit the PEM also 
exhibits a tendency to transform into another uniform or doubly modulated 
phase at low temperatures. It resembles the behaviour of the FK model in 
the strong coupling case^^The plots in Fig.[3n]are given without the effect 
of phase separation. When the segregation on the regions with different n 
values takes place, the corresponding uniform phases or both the uniform 
and the modulated ones will coexist. 



5. Tvi^o-Sublattice Pseudospin-Electron Model 

Two-sublattice 

PE y85 | 36 | 37 | 8 8176j appeared as a generalization of the usual 
PEM with the aim of more realistic description of the anharmonic subsys- 
tem of the apex oxygen ions in the YBaCaO type structures. In this case 
the pseudospin energy in the internal field is of the form hJ2i{Sii — 8^2) 
that is the reflection of th e m irror symmetry of the problem. Hamiltonian 
of the model is of the forrrPSl 



H — He + Hs + hi^s + Hs-s, 
H = He + Hs + He-s + Hs^s, 

He ^ -t^^iKi + <2) + U^{ni^nl^ + nl^ni^ 



ij sa 

Hs = -h^{S^i - 5^2) - ^^{Snl + Sn2), 

n n 



(89) 
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rr _ T \ ^ QZ QZ _ 1 \ ^ \ ^ jafi qz qz 
Hs-s — ^ 2^^nl^n2 2 2^ / . nn"-'na'-'n0 ■ 

Here, n^^ and S^^ are operators of the electron occupation number 
(S'=t,i) and pseudospin, respectively, in the n-th unit cell (a=l,2 cor- 
responds to the two apex oxygen O4 in the cell). Additionally to the term 
He-s describing the interaction between electrons and pseudospins, the 
direct interaction TJ^-s between pseudospins is included; the JS^iS^2 i^" 
teraction within the one-cell clusters is separated. Hamiltonian ([89| is in- 
variant with respect to the particle-hole transformation: nf^^ ^ 1 — n^^, 
h ^ 2g — h, ^ ^ — /i — U . 

Thermodynamics of the system described by Hamiltonian (|89p without 
a term related to the electron transfer was studied in Refs. [35l [36] . The 
mean-field approximation combined with the exact treatment of interac- 
tions within the one-cell clusters was used. For this purpose the single- 
clusters basis of states \Ri) = I^^Ji, '^fi, 'T'i2) © \Sii,Sf2), which consist 

of sixty-four state vectors, was introduced. In such a way the grand canon- 

2 

ical potential ^mf and the mean values (S^) and ^ ^ ("-no) = were 

s a—1 

calculated. 

Instead of the (S^) parameters, we can use their linear combinations: 
ri = {Sf + S2) (the order parameter for the ordered ferroelectric-like phase) 
and ^ = (5*^ — 5*1) (the parameter which is responsible for the in-phase 
reorientation of both pseudospins in the unit cell) . In the regime fi = const 
we obtain the equations for 77 and ^ from the condition of the minimum of 
the grand canonical potential $mf 

F^ — n 

(90) 



$MF = ^{(^11 + J^2)rf + (Jii - Ji2)e} (91) 
-Tin [2{el''i cosh /J i^iii^ + e-^i cosh (3 (h + iiii^^) | 

+86^^^ |e/3i cosh/j '-^^^Y'''" cosh/3f + e'^i cosh/3 (/i 

+e-^T cosh/3 (h + (^"-/i^)^ - . 
In the case n =const the condition of the minimum value of the free energy 
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FniF = ^MF + A^'T'j supplemented by the equation for chemical potential 

— « — = 

was used to determine the thermodynamically stable equilibrium states. 

When the set of equations ([90|) has a non-zero solution for i] and the 
corresponding thermodynamic potential has a minimum, then our system 
is in the pola r (ferr oelectric) phase. The obtained phase diagrams in the 
fi =const casj^^l^ are shown in Fig. [32] (all parameters are normalized 



here by Jn + J12 > 0). One can see that at 



'.^11 — J12 — 



= — 1 (the case of 



pseudospin-pseudospin interaction between different sublattices only) the 
phase transition into ferroelectric phase is of the second order. The presence 
of the intra-sublattice interaction leads to the possibility of change of the 
phase transition order and the appearance of tricritical points. An increase 
of the parameter causes the narrowing of the ferroelectric region; 

its width is also proportional to J. At J — 0, J12 = the model transforms 
into the one-sublattice PEM; the ferroelectric phase disappears and we have 
the first-order phase transition with zero value of the order parameter ry and 
a sharp change of parameter ^ = {S^ — 5*1). 



0.4 



0.0 



Tc/{Jii+Jyi 



0.3 - 



0.2 - 



0.1 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

tl/{Jll+Jl2) 




+ J12) 



-1.0 -0.5 -0.0 0.5 1.0 1.5 2.0 

V(J'll+/l2) 



Fig. 32. (Tc — h) phase diagrams at different values of parameter Jn — J12 in the regime 
= const: a) jH'^jH = -1, b) j^^" = 0. Other parameters: J/(Jii + J12) = -1, 
g/(Jll -I- J12) = 1, /i/(Jii -I- J12) = —1. The phase transitions can be of the second (solid 
lines) or of the first order (dotted lines). 



In the n=const regime, the region of h values at which the ferroelectric 
phase exists becomes broader. The corresponding phase diagrams are shown 
in Fig. [33] at different values of long range interaction J^^ — J2n' ■^nn' ^^'^ 
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fixed concentration n = 0.4. 




/l/(Jll+Jl2) ft/('^ll + '^12) 



Fig. 33. The h dependence of the temperature of ferroelectric phase transition Tc at 
different values of parameter Jn — J12 in the regime n=const: a) (Jn — Ji2)/(.^ii + 
J12) = -1, b) (Jii - Ji2)/(J'ii + J12) = 0. Other parameters: J/(Jii + J12) = 1, 
g/(Jil + J12) = 1, n = 0.4. Solid lines and dashed lines represent the second order and 
the first order phase transitions, respectively. The widely spaced dashed line corresponds 
to a one-loop approximation. Dots represent a separation area. 

For comparison, there are also shown the phase transition hnes ob- 
tained in the so-called one loop-approximation (in which the higher order 
corrections to the MFA described by two-tailed diagramJ^S! are taken into 
account). In the areas marked by points, the system is separated into two 
regions with concentrations ni and n2 (rii < n < 712). Also, one can notice 
that phase separation takes place near the border of stability region of two 
phases (the ordered phase with nonzero polarization and the disordered 
one). That is why the ordered phase spreads wider and extends up to the 
edge of the separated area. Fig. [34] illustrates such a behaviour. Here the 
dashed lines represent the region of ferroelectric type instabilities. These 
lines would separate the ferroelectric phase if there were no phase separa- 
tion. Figure [34y b) also shows that at a fixed value of asymmetry parameter 
h at concentration n < 0.75 (at low temperature), the ordered phase is 
possible only due to the phase separation. 

Now, let us focus briefly on the effect of electron transfer. As above, we 
restrict ourselves to the limit [/ — > oo. The case is considered, when the 
chemical potential is placed in the region of energy subbands (a ~ 1,2) 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



Fig. 34. Phase (T — n) diagram in a mean field approximation. The phase separation 
region is limited by solid lines. The dashed lines point to the region of ferroelectric 
instabilities. Parameters: (Jn — Ji2)/{^ii +^12) = Ij J/{Jll + Ji2) = g/{Jll+Jl2) = 1; 
a) + = 1, b) h/iJii + J12) = 1.35 . 

el\k) = sfik) = (-1)"! + + Xl% (92) 

separated by the gap equal to g. The Hubbard-I approximation in which 
the expressions (j92p are given here, was improved in Ref. j38| by the mean- 
field corrections (being of the form of the loop-like inclusions) to the elec- 
tron Green's functions. Thus, renormalization of spectrum due to the shift 
of subbands dependent on electron concentration was taken into account. 
Besides, the mean values of Hubbard operators in (|92p were determined 
self-consistently; the free energy of the pseudospin-electron systems was 
calculated in the above described GRPA approach. The pseudospin part of 
the two-sublattice Hamiltonian was taken in the MFA. 

Figure [35] illustrates the effect of electron subsystems on the shape of 
the phase diagram. With respect to the standard Mitsui model (which cor- 
responds to the g = and tij — limit), the phase boundary becomes 
asymmetric and the region of the ferroelectric phase existence is shifted to 
the higher values of h. As it has already been noted, at J12 = 0, ferroelectric 
phase does not exist and only a discontinuous change of ^ on the transition 
line takes place. Such a transition also remains in the case when Jn — 
and only the indirect interaction via electron subsystem is present. 

The changes in electron spectrum (see Refs. [37l|38]) are demonstrated 
in Fig. [36r b). In the case shown in Fig. [36l when at n =const the sep- 
aration into ferroelectric and nonpolar phases takes place, the ri{n) and 
/x(n) dependencies indicate that the ferroelectric phase appeared before the 
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separation. This is also supported by the presence of concavity in the free 
energy (dashed tangent hnes in Fig. 136( c) hnk the points with concentration 
values ni, ri2 and n^, on which the separation takes place). Hence, there 
is a separation into paraelectric and ferroelectric phases at concentrations 
ni < n < n2 and ns < n < n^. The pure ferroelectric phase exists in the 
concentration range n2 < n < n^. The separation area in the {h,n) plane 
changes its shape depending on the tjj value. In general, the electron trans- 
fer narrows the separation region''^'' The reverse effect also takes place: 
the electron spectrum is modified by a phase separation being sensitive to 
concentration value. 

One can relate the phase transitions and ferroelectric type instabilities 
described by the herein considered two-sublattice PEM to the observed di- 
electric and thermodynamic properties of the YBaCuO-type superconduct- 
ing crystals. In some experiments, YBa2Cu307_5 was found to be both 
pyroelectric and piezoelectric, implying the existence of macroscopic polar- 
ization directed along the c-axis. ■ The possibility of the existence of the 
ferroelectric-like phase was not, nevertheless, unambiguously confirmed. In 
the two- sublattice PEM the ordered polar phase is present in a rather re- 
stricted region of the model parameter values, where the ground state is 
degenerated; for example, its width along the /i-axis is determined by the 
interaction constant between the pseudospins (describing apical O4 ions) 
within the unit cell cluster. In this connection, it should be mentioned 
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Fig. 35. Dependence of temperature of the phase transition Tc on the parameter h at 
different values of interaction parameters Jn and J12 in the regime = const = —g. 
Thick lines correspond to the case tij/g = 0.2; thin lines correspond to the case tij = 0. 
a) Jii = J12 = s/2, b) Jii = g, J12 = 0. c) J\\ = J12 = 0, tij/g = 0.2. Phase transitions 
are either of the second order (solid lines) or of the first order (dashed lines). 
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Fig. 36. Dependence of the order parameter rj (a), band spectrum (b), and free energy 
(c) on electron concentration. The parameter values are: Jn = J\2 = 9/2, T/g = 0.1, 
i,j/g = 0.1, h/g = 0.5. 



that presence of the oxygen vacancies in the ch ain element of structure at 
5 > 1 effectively influences the value of field In its turn, instability 
with respect to the polar phase appearance can show up only at a certain 
nonstoichiometry (at certain values of parameter 5). 

The phase separation in the two-sublattice case is also related to the real 
structure of the YBa2Cu307_5 crystal. Besides the question of microscopic 
nature of the so-called "stripe phases" there is a problem of the genesis of 
structural inhomogeneities in a single crystal of the YBaCuO type observed 
in the experiments using Raman spectroscopjJ^ and mesoscopic structural 
investigations."" The results obtained within the model approach show the 
possibility of the Mitsui - type (due to both direct and indirect, via conduct- 
ing electrons, interactions between anharmonic structure units) mechanism 
of the mentioned instabilities. 
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6. Conclusions 

The present investigations of thermodynamics and energy spectrum of the 
PEM show a variety of phases and phase transitions. Depending on the 
thermodynamic equilibrium regimes, they manifest themselves as (i) transi- 
tions between different uniform phases, or between uniform and modulated 
phases, with the commensurate or incommensurate period of modulation 
(in the /i=const regime) or as {ii) transitions into phase separated states 
(in the n =const regime). The latter takes place when at // =const, the 
corresponding phase transition is of the first order. Such phase transitions 
can be realized at the change of temperature T, field h, chemical potential 
/i (in the case (i)) and other parameters of the model. The corresponding 
phase diagrams arc built in case of strong or weak coupling {g ^ W and 
g <W, respectively). 

A microscopic reason for phase transitions in standard PEM (without a 
direct pseudospin-pseudospin interaction) is related to the indirect effective 
coupling between pseudospins arising due to electron transfer and possess- 
ing a dynamic character. The form of such coupling depends on the electron 
concentration, temperature and the model parameter values. Consequently, 
the modulated phase appears at intermediate valuer of /i ( this corresponds 
to the electron occupation near half-filling) both at strong and at weak 
coupling. However, the formation mechanisms of the effective interaction 
are different in these cases. 

The two-sublattice PEM is a special case, close to the real HTSC systems 
of the YBaCuO type. The phase transitions described here, connected with 
the appearance of the ferroelectric state or with the jump-like change of 
the mutual orientation of pseudospins in sublattices, can have a relation to 
the segregation or bistability phenomena as well as to the development of 
ferroelectric-type instabilities in the mentioned systems. 

The analysis carried out in the weak coupling case elucidates the role 
of the transverse field 17 (having a meaning of the tunneling- like splitting 
parameter) in the obtained picture of phase transitions. As a whole, a topol- 
ogy of the phase diagrams does not change at O 7^ in comparison with 
the case = 0. The critical temperatures (including Tsc) decrease with J7. 
In the case of DOS with logarithmic singularity (at dimensionality d = 2) 
a peculiar effect is revealed: the critical temperature remains finite at any 
large values of Cl (tending to zero at O — > 00, exclusively). This is the dif- 
ference with respect to the behaviour of the systems with direct interaction 
(e.g., Ising model with transverse field). The important property is that at 
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^ the superconducting phase can appear in the PEM. Such a phase 
competes with the modulated one and is stable at the electron occupancy 
near the upper (or lower) edge of the electron band. 

As is seen from the results obtained in DMFT for the simplified model, 
the structure of electron spectrum of the PEM is different in cases g ^ W 
and g < W: the split subbands due to interaction (even at J7 = 0) or a single 
band, respectively. The similar spectrum is obtained in the approximations, 
based on the Hubbard-I scheme (e.g. GRPA) in the first case or the Hartree- 
Fock ap proach in the second one. As was shown by Zeyher, Kulic, and 
Gehlhoffj I I I GRPA keeps in a systematic way all terms of the leading 
order of a so-called expansion, where N is the local spin component 
number on a lattice site. Though such approximations do not adequately 
reproduce all the features of spectrum, the obtained phase diagrams are 
in a good agreement with the ones constructed in DMFT. An interesting 
feature of spectrum is that there exists a critical value of g: at g > gc & 
gap appears (in the case of simplified PEM) and a metal-insulator type 
transition takes place. 

Pseudospin-electron model (PEM) can be considered as a generalization 
of the Falicov-Kimball (FK) model to the case of different thermodynamic 
equilibrium regimes as well as an extension of the latter model due to the 
inclusion of the pseudospin dynamics and the Hubbard type correlations. 
The PEM possesses a similar variety of phase transitions but there are 
differences in the conditions of their realization and in the criteria of the 
appearance of different phases. In the above considered cases such differ- 
ences are discussed and a comparison with the results for the FK model is 
made. 

From the point of view of theoretical studies, the investigations of the 
PEM are far from complete. Another interesting problem is connected with 
the thermodynamics of the PEM with electron transfer dX U ^ Q and 
^ 0. The investigations performed revealed only the existence of insta- 
bilities connected with certain values of the wave vector q but the phase 
diagrams determining the regions of existence of different phases have not 
been built so far. The intermediate coupling case {g ~ W), that was not 
adequately investigated even for the simplified PEM, calls for more detailed 
consideration. An important point is to complete the study of the collec- 
tive excitation spectrum (connected with pseudospin reorientation, electron 
transitions and polaron effect) and the dynamic susceptibility. Among the 
possible generalizations of the model, one can note an extension to the 
cases with the asymmetric electron transfer (in the spirit of the asymmet- 
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ric Hubbard model) and with pseudospin S > 1/2. They appear to be quite 
promising in connection with investigations of the ionic transport based on 
the lattice models as well as in the study of the ion intercalation processes 
(see e.g. Ref . [84l for a recent development in this field). 

Acknowledgements 

The author is grateful to A. Shvaika and T. Mysakovych for reading the 
manuscript and technical help. 

References 

1. K. A. Miiller, Z. Phys. B 80, 193 (1990). 

2. S. Conradson and I. D. Raistrick, Science 243, 1340 (1989). 

3. J. Mustre de Leon, S. D. Conradson, L. Batistic, A. R. Bishop, I. D. 
Raistrick, M. C. Aronson, and F. H. Garzon, Phys. Rev. B 45, 2447 (1992). 

4. D. Mihajlovic and C. M. Foster, Solid State Commun. 74, 753 (1990). 

5. G. Ruani, C. Taliami, M. Muccini, K. Conder, E. Kaldis, H. Keller, D. Zech, 
and K. A. Miiller, Physica C 226, 101 (1994). 

6. M. N. Iliev, V. G. Hadjiev, and V. G. Ivanov, Journ. Raman Spectr. 27, 
333 (1996). 

7. N. Poulakis, D. Palles, E. Liarokapis, K. Conder, E. Kaldis, and K. A. 
Miiller, Phys. Rev. B 53, R534 (1996). 

8. L. R. Testardi, W. G. Moulton, H. Matials, H. K. Ng, and C. M. Rey, Phys. 
Rev. B 37, 2324 (1988). 

9. V. Miiller, C. Hucho, and D. Maurer, Ferroelectrics 130, 45 (1992). 

10. J. Mustre de Leon, L Batistic, A. R. Bishop, S. D. Conradson, and S. A. 
Trugman, Phys. Rev. Lett. 68, 3236 (1992). 

11. J. Ranninger and U. Thibblin, Phys. Rev. B 45, 7730 (1992). 

12. A. P. Saiko and V. E. Gusakov, JETP 108, 757 (1995). 

13. P. Schweiss, W. Reichardt, M. Braden, G. Collin, G. Heger, H. Glaus, and 
A. Erb, Phys. Rev. B 49, 1387 (1994). 

14. C. Ambrosch-Draxl and P. Knoll Physica B 194-196, 2091 (1994). 

15. J. Rohler, in: Materials and Crystallographic Aspects of HTc-su- 
perconductivity, edited by E. Kaldis (Kluwer, Dordrecht, 1994), p. 353. 

16. R. J. Cava and A. W. Hewat, Physica C C165, 419 (1990). 

17. J. R. Hardy and J. W. Flocken, Phys. Rev. Lett. 60, 2191 (1988). 

18. A. Bussman-Holder, A. Simon, and H. Biittner, Phys. Rev. B 39, 207 (1989). 

19. N. M. Plakida, Physica Scripta 29, 77 (1989). 

20. J. E. Hirsch and S. Tang, Phys. Rev. B 40, 2179 (1989). 

21. N. M. Plakida and V. S. Udovenko, Mod. Phys. Lett. B6, 541 (1992). 

22. M. Frick, W. van der Linden, I. Morgenstern, and H. Raedt, Z. Phys. B 81, 
327 (1990). 

23. I. V. Stasyuk, A. M. Shvaika, and E. Schachinger, Physica C 213, 57 (1993). 

24. I. V. Stasyuk and A. M. Shvaika, Fiz. Nizk. Temp. 22, 535 (1996). 



February 2, 2008 3:3 



WSPC/Trim Size: 9in x 6in for Review Volume 



stasyuk'new 



Phase Transitions in the Pseudospin- Electron Model 59 

25. I. V. Stasyuk and A. M. Shvaika, Ferroelectrics 192, 1 (1997). 

26. I. V. Stasyuk and A. M. Shvaika, Condens. Matter Phys. 3, 134 (1994). 

27. V. Koorting, Q. Yuan, P.J. Hirschfeld, T. Kopp, J. Mannhart, Phys. Rev. 
B 71, 104510 (2005). 

28. N. Pavlenko, T. Kopp, Phys. Rev. B 72, 174516 (2005). 

29. I. V. Stasyuk and A. M. Shvaika, Journ. Phys. Studies, 3, 177 (1999). 

30. I. V. Stasyuk, A. M. Shvaika, and K. V. Tabunshchyk, Condens. Matter 
Phys. 2, 109 (1999). 

31. I. V. Stasyuk, A. M. Shvaika, and K. V. Tabunshchyk, Ukrainian Journ. of 
Phys. 45, 520 (2000). 

32. I. V. Stasyuk and T. S. Mysakovych, Journ. Phys. Studies 5, 268 (2001). 

33. I. V. Stasyuk and T. S. Mysakovych, Condens. Matter Phys. 5, 473 (2002). 

34. T. S. Mysakovych and I. V. Stasyuk, Ukrainian Journ. of Phys. 49, 607 
(2004). 

35. O. D. Danyliv and I. V. Stasyuk, Condens. Matter Phys. 7, 163 (1996). 

36. O. D. Danyhv, Physica C 309, 303 (1998). 

37. O. D. Danyliv and I. V. Stasyuk, Condens. Matter Phys. 5, 523 (2002). 

38. I. V. Stasyuk and O. D. Danyliv, Phys. Stat. Sol. 219, 229 (2000). 

39. Yu. A. Izyumov and B. M. Letfulov, J. Phys.: Cond. Matter 2, 8905 (1990). 

40. I. V. Stasyuk and T. S. Mysakovych, Journ. Phys. Studies 3, 344 (1999). 

41. I. V. Stasyuk and T. S. Mysakovych, Physica C 341-348, 171 (2000). 

42. T. S. Mysakovych and I. V. Stasyuk, Condons. Matter Phys. 7, 347 (2004). 

43. J. K. Freericks and V. Zlatic, Rev. Mod. Phys. 75, 1333 (2003). 

44. I. V. Stasyuk, R. Ya. Stetsiv, and Yu. V. Sizonenko, Cond. Matter. Phys. 
5, 685 (2002). 

45. I. V. Stasyuk and Yu. I. Dublenych, Phys. Rev. B 72, 224209 (2005). 

46. W. Wagemaker, G. J. Kearby. A. A. van Well, H. Mutka, and F. M. Mulder, 
J. Am. Chem. Soc. 125, 840 (2003). 

47. V. Koerting, Qingshan Yuan, P. J. Hirschfeld, T. Kopp, and J. Mannhart, 
Phys. Rev. B 71, 104510 (2005). 

48. W. Metzner and D. VoUhardt, Phys. Rev. Lett. 62, 324 (1989). 

49. A. Georges, G. Kotliar, W. Krauth, and M. J. Rosenberg, Rev. Mod. Phys. 
68, 13 (1996). 

50. W. Metzner, Phys. Rev. B 43, 8549 (1991). 

51. U. Brandt and C. Mielsch, Z. Phys. B 75, 365 (1989); iMd. 79, 295 (1990). 

52. G. Baym and L. P. Kadanoff, Phys. Rev. 124, 287 (1961); G. Baym, Phys. 
Rev. 127, 1391 (1962). 

53. E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

54. V. J. Emery and S. A. Kivclson, Physica C 209, 597 (1993); U. Low, V. J. 
Emery, K. Fabrieius, and S.A. Kivelson, Phys. Rev. Lett. 72, 1918 (1994). 

55. T. Mertelj, J. Demsar, B. Podobnik, L Poberaj, and D. Mihailovic, Phys. 
Rev. B 55, 6061 (1997). 

56. L V. Stasyuk and A. M. Shvaika, Czech. J. Phys. 46, 961 (1996). 

57. J. K. Freericks, Phys. Rev. B 47, 9263 (1993). 

58. P. M. Slobodyan and I. V. Stasyuk, Theor. Math. Phys. USSR 19, 616 
(1974); ( Teor. Mat. Fiz. 19, 423 (1974)). 



February 2, 2008 3:3 WSPC/Trim Size: 9in x 6in for Review Volume stasyuk'new 



60 /. Stasyuk 

59. Yu. A. Izyumov, F. A. Kassan-Ogly and Yu. N. Skryabin, Field Methods in 
the Theory of Ferromagnetism (Nauka, Moscow, 1974) (in Russian) . 

60. I. V. Stasyuk and Yu. Havrylyuk, Condons. Matter Phys. 2, 487 (1999). 

61. J. K. Freericks and R. Lcmanski, Phys. Rev. B 61, 13438 (2000). 

62. J. K. Freericks, Ch. Gruber, and N. Macris, Phys. Rev. B 60, 1617 (1999). 

63. U. Brandt and C. Mielsch, Z. Phys. B 82, 37 (1991). 

64. B. M. Letfulov, Eur. Phys. J. 4, 447 (1998); 11, 423 (1999). 

65. I. V. Stasyuk and A. M. Shvaika, Acta Phys. Polon. A 84, 293 (1993). 

66. J. K. Freericks, Phys. Rev. B 48, 14797 (1993). 

67. N. M. Plakida and V.S. Udovenko, Sverkhpr. Fiz. Khim. Tekhn. 5, 775 
(1992). 

68. Yu. A. Izyumov, M. I. Katsnel'son and Yu. N. Skryabin, Magnetism of 
Itinerant Electrons (Nauka, Moscow, 1992) (in Russian). 

69. J. K. Freericks, M. Jarrell, and D. J. Scalapino Phys. Rev. B 48, 6302 (1993); 
ibid. B50, 403 (1994). 

70. C. S. Owen and D. ,J. Scalapino, Physica 55, 691 (1971). 

71. W. P. Su, Phys. Rev. B 67, 092502 (2003). 

72. N. N. Bogolyubov, V. L. Aksenov, and N. M. Plakida, Teor. Mat. Fiz. 93, 
371 (1992) (in Russian). 

73. I. V. Stasyuk and A. M. Shvaika, Acta Phys. Polon. A 85, 363 (1994). 

74. I. V. Stasyuk, A. M. Shvaika, and E. Schachinger, Physica B 194-196, 1965 

(1994) . 

75. I. V. Stasyuk and A. M. Shvaika, Physica C 235-240, 2173 (1994). 

76. I. V. Stasyuk, A. M. Shvaika, and O. D. Danyliv, Molec. Phys. Rep. 9, 61 

(1995) . 

77. X. X. Xi, C. Doughty, A. Walkcnhorst, C. Kwon, Q. Li, and T. Venkatesan, 
Phys. Rev. Lett. 68, 1240 (1992). 

78. D. Mihailovic and A. J. Heeger, Solid State Comm. 75, 319 (1990). 

79. I. V. Stasyuk and O. V. Velychko, Ukrainian Journ. of Phys. 44, 772 (1999). 

80. V. M. Browning, Phys. Rev. B 65, 2860 (1997). 

81. M. L. Kulic, R. Zcyhcr, Phys. Rev. B 49, 4395 (1994). 

82. L. Gehlhoff, R. Zeyher, Phys. Rev. B 52, 4635 (1995). 

83. R. Zeyher, M. L. Kulic, Phys. Rev. B 53, 2850 (1996). 

84. F. Zhou, T. Maxisch, G. Ceder, Phys. Rev. Lett. 97, 155704 (2006). 



2, 2008 3:3 



WSPC/Trim Size: 9in x 6in for Review Volume 



stasyuk': 



INDEX 



approximation 

coherent potential (CPA), 6 
generalized random phase (GRPA), 
12 

Hubbard-I, 8 
mean-field, 26 
Migdal-Ehashberg, 39 

charge density wave (CDW), 19 
correlation 
electron U, 2 

density of states (DOS), 6 
dielectric 

instability, 47 

response, 47 

susceptibility, 45 

equation 

Bethe-Salpeter, 39 

function 
Green's, 5 

Matsubara, 35 

Zubarev, 35 
pair correlation, 15, 34 

Hamiltonian 
Hubbard, 2 

model 

Falicov-Kimball (FK), 4 
Hubbard, 2 
Mitsui, 53 

pseudospin-electron (PEM), 2 



simplified, 12 
two-sublattice, 48 

phase 

chess-board, 20 
coexistence, 27 
diagram, 8, 34 
ferroelectric, 50 
incommensurate, 36 
separation, 10 
transition 

first-order, 8 
Mott, 7 

point 

spinodal, 36 

tricritical, 32, 50 
pseudospin, 2 

response 

"isolated", 35 
isothermal, 34 

spinodal, 36 

state 

superconducting (SC), 38 
theory 

dynamical mean field (DMFT), 5 
tunneling splitting, 2, 27, 41 

YBaCuO-type superconducting 
crystals, 2, 35 



61 



